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The  subject  of  this  annual  report  is  constitutive  modeling  of 
cohesionless  soil,  for  both  laboratory  standard  static  test  conditions  and 
insitu  impulsive  dynamic  load  conditions.  Laboratory  test  results  are 
used  to  quantify  the  parameters  for  several  different  types  of 
constitutive  models.  The  ability  of  these  models  to  predict  laboratory 
behavior  Is  then  compared  and  evaluated.  The  modeling  of  insitu  soil 
response  to  explosive  events  is  considered,  and  the  laboratory-derived 
models  are  tested  for  their  convenience  and  accuracy  in  predicting  ground 
motions.  Several  phenomena  which  occur  in  the  laboratory  and/or  insitu, 
and  which  may  be  Important  in  certain  cases  but  were  not  reproduced  by  the 
models  used,  are  discussed.  An  "ideal"  testing  device  is  proposed,  which 
can  exercise  the  Important  stress-strain  mechanisms  of  explosively  loaded 
cohesionless  soils.  The  complexity  of  a  constitutive  model  which  Includes 
all  the  above  mechanisms  is  considered. 

The  following  Is  the  Statement  of  Work  for  this  effort,  as  it  appears 
In  Contract  Document  No.  F49620-80-C-0088,  FY82  Modification: 

a.  Use  the  Soil  Element  Model  to  study  the  response  of  soils  to 
laboratory  boundary  conditions,  utilizing  a  variety  of  dynamic 
soil  models.  Emphasis  will  be  placed  on  understanding  each 
model's  capability  to  reproduce  a  large  suite  of  laboratory  test 
data  and  more  complex  stress/ strain  paths  representing  insitu 
dynamic  loadings. 

b.  Using  the  results  from  (a),  characterize  the  "ideal"  laboratory 
test(s)  required  to  describe  insitu  response  of  soils  to  complex 
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dynamic  loadings.  Compare  the  resulting  stress/strain  paths  with 
those  of  current  laboratory  tests. 

c.  Continue  the  CISTFTON  study  of  Insltu  soil  response.  Begin 
efforts  to  model  the  response  of  at  least  one  soil  for  which 
Insltu  testing  and  large  explosive  event  data  are  available.  The 
purpose  Is  to  study  the  process  of  laboratory  tests,  insitu 
tests,  model  fitting  and  calculation  of  a  dynamic  event  to 
confirm  deficiencies  In  current  models. 

d.  Begin  efforts,  based  on  a,  b,  c,  to  outline  a  new  approach  for 
constitutive  modeling  for  soils. 


j 
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2.0  PROGRESS  SUMMARY 

2.1  Soil  Element  Model  (SEM)  Development 

The  SEM  Isa  FORTRAN  computer  program  which  was  written  during  the 
Jast  reporting  period  (Ref.  1)  to  examine  some  of  the  questions  regarding 
dynamic  soil  modeling.  It  consists  of  a  driver,  a  set  of  boundary 
condition  algorithms,  and  a  set  of  constitutive  relationships.  The  SEM 
allows  comparison  of  material  models  when  they  are  subjected  to  various 
kinds  of  static  and  dynamic  stress  and/or  strain  paths.  It  can  also  be 
used  for  parameter  variation  studies,  and  as  a  model  parameter  fitting 
tool. 

During  this  reporting  period,  the  SEM  continued  to  be  the  primary 
tool  for  constitutive  model  evaluation  and  development.  Two  additional 
models  were  added  to  those  already  being  studied:  Duncan  and  Chang's 
hyperbolic  curve-fitting  model  (Ref.  2),  and  a  standard  linear 
viscoelastic  model.  Also,  Lade's  model  for  cohesionless  soil  (Ref.  3)  has 
been  programmed  and  is  currently  being  tested. 

Several  models  already  installed  were  refined  to  improve  their 
behavior  under  various  loading  conditions.  The  AFWL  Engineering  model 
hydrostat  was  changed  to  a  pressure  versus  volumetric  strain  formulation 
In  place  of  a  pressure  versus  excess  compression  relationship.  This 
resolved  some  inconsistencies  between  assumed  incrementally  elastic 
properties.  A  multilinear  unload-reload  relation  was  added  to  the  AFWL 
Engineering  model,  which  increased  the  amount  of  strain  recovery  and 
allowed  for  hysteresis  loops.  The  AFWL  Engineering  model  failure  surface 
was  also  modified  to  give  different  shear  strengths  in  triaxial 
compression  and  extension,  to  produce  a  Mohr-Coulomb-like  cross-section  in 
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the  octahedral  plane.  Tensile  behavior  In  several  models  was  studied  and 
In  most  cases  improved.  It  was  found  that  in  many  loading  situations  the 
response  to  tensile  stress  significantly  affected  overall  behavior. 

Various  forms  of  tensile  stress  cutoff  and  volumetric  strain  tracking  have 
been  tested  In  the  SEM. 

The  versatility  of  the  SEM  in  subjecting  a  particular  model  to 
various  stress/straln  paths  was  greatly  Improved  during  the  past  year  by 
addition  of  a  one-dimensional  wave  propagation  feature.  The  finite 
difference  logic  for  this  feature  was  adapted  from  SNEAKY  (Ref.  4),  and 
can  use  any  of  the  SEM  stress- strain  models.  The  code  can  treat  planar, 
cylindrical,  and  spherical  geometries,  and  accommodates  up  to  ten 
different  material  layers,  each  of  which  can  have  a  different  constitutive 
model.  Although  restricted  to  one  dimension.  Inclusion  of  wave 
propagation  in  the  SEM  is  an  Important  step  because  it  makes  model 
evaluation  under  Impulsive  loading  much  easier,  and  the  results  directly 
applicable  to  finite  difference  techniques. 

Improved  plotting  has  made  assessment  of  results  from  SEM  studies 
easier  and  quicker.  The  curves  resulting  from  several  problems  can  be 
overlaid  for  qualitative  comparison.  Quantitative  comparisons  can  be  made 
by  calculating  the  area  between  two  curves  as  a  continuous  residual.  In 
addition,  any  quantity  resulting  from  a  wave  propagation  calculation  can 
be  plotted  versus  time,  another  time  varying  quantity,  or  grid  position. 

2.2  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  Testing 
Program 

Laboratory  tests  have  been  conducted  on  two  recompacted  sands,  Misers 
Bluff  sand  and  Reid  Bedford  Model  sand.  The  data  already  available  for 
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these  materials,  plus  the  above  test  results  comprise  a  fairly  complete 
and  accurate  definition  of  both  sands'  uniaxial  strain  and  standard 
trl axial  behavior  under  static  and  dynamic  (rapid)  loading.  WES  has 
published  three  reports  (Refs.  5,  6,  and  7)  which  document  the  testing 
^program  and  provide  Initial  analyses  of  the  data. 

2.3  Application  of  Laboratory  Data  In  Modeling 

The  data  from  the  WES  program  was  utilized  to  study  and  model  the 
behavior  of  coheslonless  sands.  The  SEM  was  used  to  arrive  at  model  fits 
to  the  data  for  most  of  the  Implemented  constitutive  relations.  The 
ability  of  each  model  to  predict  laboratory  responses  was  evaluated  and 
compared.  This  work  Is  described  In  detail  In  Section  4.0  of  this  report. 

2.4  Strain  Path  Response 

Model  response  comparisons  were  made  under  the  "typical"  strain  path 
excursions  postulated  by  Workman,  et  al.  (Ref.  8).  These  strain  paths  are 
those  expected  from  a  truly  spherical  field  of  motion  resulting  from  a 
burled  explosion.  The  primary  purpose  of  this  exercise  was  to  see  If  the 
models  Installed  In  the  SEM  predicted  reasonable  stress  paths  when 
exercised  over  the  postulated  strain  paths.  It  was  found  that  each  model 
produced  a  reasonable  stress  path,  and  that  the  most  important  material 
parameters  affecting  the  results  were  shear  strength  and  tensile  behavior. 

A  complete  study  of  the  strain  path  response  technique  for  modeling 
material  response  has  not  been  made.  However,  some  Initial  observations 
can  be  made.  The  method  requires  a  stress  path  resulting  from  each  strain 
path.  Aside  from  the  testing  difficulties  Involved  In  measuring  such 
response  In  the  laboratory,  the  representative  strain  paths  may  not  always 
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be  unique,  especially  under  conditions  of  tensile  or  shear  failure.  Under 
these  circumstances  It  would  be  difficult  to  utilize  the  resulting 
stress-strain  relations  as  a  complete  material  model. 

2.5  Wave  Propagation  Studies 

The  models  derived  from  the  laboratory  results  were  used  to  predict 
ground  motions  due  to  blast  loading  for  uniaxial  strain,  plane-strain,  and 
axl symmetrl c  and  spherically  symmetric  geometries.  The  wave  propagation 
studies  revealed  much  about  the  models  which  was  not  Immediately  obvious 
from  laboratory  behavior  predictions,  and  forced  them  to  become  completely 
functional  In  a  finite  difference  code.  The  dependency  of  model  response 
on  calculatlonal  details  such  as  time  step,  grid  size,  and  artificial 
viscosity  was  explored.  For  example,  it  was  quickly  confirmed  that  the 
cap  model  may  require  subcycling  within  a  Courant-conditlon-defined  time 
step  to  ensure  correct  model  response.  It  was  also  found  that  the 
spherical  and  cylindrical  one-dimensional  wave  propagation  calculations 
exercised  several  Important  aspects  of  the  models  more  completely  than  did 
the  planar  calculations.  Expansive  strains  are  produced  In  the 
cylindrical  and  spherical  calculations  which  test  the  tensile 
relationships,  and  large  stress  differences  are  generated  which  test  the 
adequacy  of  shearing  behavior. 

Calculated  model  responses  were  compared  with  available  measured 
Insltu  responses.  In  some  cases,  the  two  did  not  compare  well  and  model 
parameters  were  adjusted  to  give  an  insltu-based  model.  Observing  the 
differences  between  laboratory-based  and  insitu-based  models  can  help 
discern  errors  In  material  parameters  based  on  a  laboratory  testing 
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2.6  Interactions 

Appendix  D  contains  a  paper  presented  at  the  International  Conference 
on  Constitutive  Laws  for  Engineering  Materials,  In  Tuscon,  Arizona  on 
13  January  1983,  and  published  In  the  Proceedings.  The  paper  discusses 
.the  goals  of  this  research  and  briefly  summarizes  some  of  the  work 
described  above. 

Also  Included  In  Appendix  D  Is  the  abstract  of  a  paper  accepted  for 
presentation  at  the  Symposium  on  Interaction  of  Non-Nuclear  Munitions  with 
Structures,  to  be  held  at  the  USAF  Academy  from  9-13  May  1983.  That  paper 
describes  the  direct  application  of  methods  developed  under  this  effort  to 
an  actual  site  and  laboratory/lnsltu  test  program. 
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3.0  LABORATORY  BEHAVIOR  OF  SAND 

3.1  Introduction 

Soil  behavior  observed  In  laboratory  tests  Is  generally  limited  in 
Its  nature  by  several  factors.  First,  laboratory  response  is  physically 
constrained  by  test  Item  size  and  testing  apparatus  boundaries.  Second, 

T 

the  variety  of  stress  (or  strain)  paths  to  which  a  test  Item  can  be 
subjected  Is  limited  by  test  apparatus  capability.  Third,  principal 
stress  and  strain  directions  are  constant,  and  coincide  with  the  test 
specimen  axes  in  the  uniaxial  and  trlaxlal  apparatuses,  which  are  the 
mainstay  of  standard  soil  testing.  In  addition,  due  to  the  cylindrical 
nature  of  the  samples,  two  of  the  principal  stresses  and  strains  are 
always  equal.  Nevertheless,  practical  considerations  dictate  continued 
reliance  on  standard  laboratory  tests  to  develop  constitutive 
relationships  for  most  geologic  material s. Many  assumptions  are  then 
required  to  generalize  the  test  results  to  predict  response  under  more 
general  states  of  stress  and  strain. 

3.2  Material  Description 

The  dynamic  stress-strain  behavior  of  three  cohesionless  sands  was 
studied,  using  constitutive  models  based  on  laboratory  test  data.  The 
three  sands  were:  Misers  Bluff  (MB)  sand,  Reid  Bedford  Model  (RB)  sand, 
and  Ralston  Valley  (RV)  alluvium.  The  MB  and  RB  sands  were  tested  In 
conjunction  with  this  project,  and  will  receive  the  most  attention.  RV 
alluvium  results  are  from  a  separate  testing  effort  (Refs.  9  and  10),  and 
are  of  Interest  because  of  the  large  amount  of  complementary  insltu  data 
also  available  for  this  material. 
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MB  sand  Is  a  medium  to  coarse  grained  silty  sand  obtained  at  Planet 
Ranch,  Arizona.  Planet  Ranch  was  the  site  of  the  Misers  Bluff,  Phase  II 
HE  experiments  in  1978,  and  the  CIST  19  material  properties  test  In  1977. 
RB  sand  Is  a  fine  grained,  uniform  sand  obtained  along  the  Big  Black  River 
In  Mississippi.  It  Is  often  used  by  WES  for  test  apparatus  development 
and  calibration  work  because  of  Its  easy  availability,  uniformity,  and 
fairly  well  defined  properties.  RV  alluvium  is  a  medium  to  coarse  grained 
silty  sand  obtained  from  the  top  six  meters  of  the  Ralston  Valley 
geotechnical  study  site  between  Tonopah,  Nevada  and  the  northwest  corner 
of  the  Nevada  Test  Site.  Average  grain  size  distribution  curves  for  the 
three  sands  are  shown  In  Figure  1.  Average  physical  properties  are  given 
In  Table  1. 

3.3  Laboratory  Tests  Performed 

Existing  data  was  first  collected  for  dry  RB  sand  and  saturated  MB 
sand  (Ref.  5).  Uniaxial  compression  and  trlaxial  shear  tests  were  then 
performed  on  saturated  RB  sand  and  dry  MB  sand  to  complete  the  data  base. 

A  brief  description  of  each  laboratory  test  performed  at  WES  Is  given  In 
Appendix  A  (Refs.  5  and  6).  Enough  tests  were  done  so  that  representative 
relations  could  be  determined  by  weighting  and  averaging  all  test  results 
(Fig.  2).  However,  there  was  not  a  sufficient  number  of  tests  done  to 
warrant  a  statistical  determination  of  mean  curves  or  standard 
deviations.  Table  2  shows  the  number  of  each  type  of  test  done  on  MB  and 
RB  sand  during  both  previous  WES  efforts  and  the  current  ARA/WES  effort. 

The  laboratory  response  of  RV  sand  Is  less  well  defined  than  for  MB 
or  RB  sand,  but  several  static  and  (tynamlc  uniaxial  compression  and 
trlaxial  shear  tests  were  done  on  samples  from  the  depth  range  of 
Interest.  Representative  curves  have  been  estimated  from  these  results. 
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Figure  1.  Average  Grain  Size  Distribution  Curves  for  MB,  RB,  and  RV  Sands 


Table  1.  Average  Physical  Properties  of  MB,  RB,  and  RV  Sands 


MB  Sand  Silty  Sand  (SW-SM)  Brown 


2 

Volume  (cm) 


RB  Sand  Sand  (SP)  Brown 


Volume  (cm3) 


Gs  =  2,67 
e  =  0.57 

w  =  21% 

S  =  100% 

Yt  -  2.07 

w  =  0.4% 

..  n  S  =  1.9% 

Air  Dry: 

Yd  =  1.72  g/cc 


G$  =  2.65 
e  =  0.61 


w  =  23% 

S  =  100% 

Yt  s  2.03  g/cc 

w  =  0.1% 


Air  Dry: 


S  =  0.3% 

Yd  =  1.65  g/cc 


RV  Alluvium  Silty  Sand  with  Gravel  (SW-SM)  Brown  q  =2.63 


Figure  2.  Example  Determination  of  a  Representative  UX  Compression  Curve 


Table  2.  Laboratory  Tests  on  MB  and  RB  Sands 


Saturated 

Drained  Undrained 


Isotropic  Compression 
(Static) 


Uniaxial 

Strain 


Static 


Dynamic 


Uniaxial  Stra1n/K0  4 


Uniaxial  Straln/Null  0 


Trl axial  Compression  6 


Trlaxlal  Extension  3 


6*  2* 


23*  5* 


♦Tests  performed  under  an  earlier  WES  Effort. 

♦See  Appendix  A  for  descriptions  of  test  procedures. 


Saturated 

Drained  Undrained 


Three  types  of  tests  were  conducted  on  remolded  MB  and  RB  sand: 
undralned-dry  (UD),  dralned-saturated  (DS),  and  undralned-saturated  (US). 
The  Initial  void  ratio  was  the  same  for  each  sand.  Results  discussed 
below  are  only  a  summary  In  the  form  of  representative  curves.  References 
6  and  7  should  be  consulted  for  the  complete  data  set. 

3.4  Laboratory  Test  Results 

Uniaxial  compressibility  Is  shown  In  Figure  3  for  MB  and  RB  sand. 
These  representative  curves  combine  both  static  and  dynamic  test  results. 
The  curves  marked  UD  are  not  undrained  In  the  conventional  sense  ( 1 . e. , 
saturated  with  constant  water  content).  If  that  were  the  case  there  would 
be  little  deformation  In  a  uniaxial  compression  test.  As  Appendix  A 
Indicates,  a  UD  test  is  a  test  on  an  air  dry  specimen  with  pore  air  flow 
prevented.  The  apparent  compressibility  difference  between  the 
representative  UD  and  DS  curves  may  be  due  to  the  same  water  lubricating 
effect  which  also  Influences  compaction  test  results.  However,  it  must  be 
noted  that  the  representative  UD  and  DS  curves  were  each  constructed  by  a 
Judgmental  process,  from  sets  of  curves  exhibiting  far  greater 
compressibility  differences  than  that  between  the  two  resulting 
representative  curves.  Thus,  the  apparent  compressibility  difference 
between  the  representative  UD  and  DS  curves  may  not  be  significant.  WES 
reported  the  difference  because  It  was  observed.  The  difference  is 
acknowledged  here  because  associated  with  it  are  two  types  of  Initial 
compressibility  relations  (concave  and  convex  to  the  strain  axis)  which 
both  need  to  be  constltutlvely  modeled. 

MB  and  RB  sands  are  nonlinear,  fairly  compressible  In  the  UD  or  DS 
conditions,  and  recover  only  a  relatively  small  amount  of  strain  upon 
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unloading.  At  lower  stress  levels  (below  10  MPa)  the  DS  response  Is 
apparently  stlffer  than  the  UD  response.  Above  10  MPa  both  the  UD  and  DS 
samples  show  similar  response.  The  US  uniaxial  strain  response  is  linear 
and  elastic,  with  a  modulus  predictable  from  mixture  theory,  using  the 
bulk  modulus  of  water,  soil  porosity,  and  an  assumed  Poisson's  Ratio  for 

T 

the  water-sand  mixture.  For  MB  sand,  with  n  =  0.36  and  with  the  bulk 
modulus  of  water,  «w  =  2070  MPa,  the  constrained  modulus  of  the 
water- sand  mixture,  M$,  may  be  computed  as  follows: 


M 


s 


3(1  -  v)  Kw 
TT+T)"  “n 


Ms  = 


3(1  - 

— nr 


.49)2070 

♦  .47} 


OF 


=  5900  MPa 


This  compares  well  with  the  observed  representative  value  of  6000  MPa. 

Also  worth  noting  Is  the  UD  and  DS  hysteretlc  unload-reload  behavior.  The 
observed  trend  is  for  these  loops  to  become  larger  when  unloading  and 
reloading  from  higher  stress  levels.  This  effect  is  particularly 
Important  when  considering  cyclic  loading,  bt.ause  of  the  energy 
dissipation  it  causes. 

Figure  4  is  a  representative  UX  compressibility  curve  for  RV  alluvium 
(Ref.  11),  based  on  drained  tests  on  undisturbed  samples  at  their  natural 
water  content.  The  compressibility  of  RV  sand  is  very  similar  to  that  of 
MB  sand. 

UX  stress  paths  from  UX/Kq  tests  are  shown  for  MB  and  RB  sand  in 
Figure  5.  For  both  materials  the  US  responses  were  very  close,  as  were 
the  DS  responses.  Based  on  elastic  theory,  the  slope  of  the  UX  stress 
path,  s  *  a(oj-o3)/ap,  Is  directly  related  to  Poisson's  Ratio,  v,  by 

the  relation: 
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Vertical  Stress  (MPa) 


V 


(3.1) 


3(Aoj  -  A03)  3(1  - 

AOj  +  “  j  + 


-v 

77 


*  3(J 


T+" 


2v, 


This  equation  Is  plotted  In  Figure  6A.  Figure  6B  shows  the  observed 
variation  In  Poisson's  Ratio  for  MB  and  RB  sands  as  calculated  from  the 
ratio  of  radial  to  axial  strain  Increments.  Poisson's  Ratio  Increases 
from  about  0.25  to  0.45  on  loading,  displays  a  sharp  discontinuity  at  the 
start  of  unloading,  and  then  increases  again  from  0.15  to  0.25.  For  the 
US  case,  the  water-sand  mixture  Is  very  nearly  Incompressible,  and  loads 
and  unloads  at  v  »  0.49. 

The  standard  tri axial  shear  response  of  MB  and  RB  sands  was  evaluated 
for  three  constant  confining  pressures:  0.14,  1.72,  and  3.45  MPa.  A  few 
additional  tests  were  done  on  RB  sand  at  a  constant  confining  pressure  of 
6.90  MPa.  As  In  the  uniaxial  tests,  the  sand  was  tested  under  three 
conditions:  undralned-dry  (UD),  drained-saturated  (DS),  and  undrained- 
saturated  (US).  The  representative  UD  and  DS  tri axial  responses  for  MB 
sand  were  very  similar,  and  In  Figure  7  the  two  curves  are  overlayed.  The, 
US  responses  are  also  shown  In  Figure  7.  When  the  DS  and  US  results  are 
compared.  It  Is  evident  that  effective  stress  controls  behavior  In  the 
undrained  state.  The  slope  of  the  US  effective  stress  path  indicates 
whether  the  sample  Is  attempting  to  dilate  (and  therefore  generate 
negative  excess  pore  pressure)  or  compress  (and  therefore  generate 
positive  excess  pore  pressure).  A  positive  slope  less  than  3  on  1 
Indicates  a  tendency  to  dilate;  a  positive  slope  greater  than  3  on  1  or  a 
negative  slope  Indicates  a  tendency  to  compress.  This  Is  because  when 
o3  Is  constant,  as  It  was  In  the  above  tests,  then: 
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Poissons  Ratio 


S,  Slope  of  UX  Stress  Path 
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Representative  TX  Stress-Strain  and  Effective  Stress  Path  Relations  for  MB  Sand 


3 


AOj  -  Aoj 


AOj  +  2 Aoj 


Aaj  -  3&u 


Clearly  the  coupling  between  shear  strain  and  volumetric  strain  Is  a 
fundamental  soil  stress-strain  mechanism  which  must  be  constltutively 

* 

modeled. 

In  order  to  explore  the  symmetry  of  the  failure  surface  about  the 
p-axls  In  the  trlaxlal  plane,  three  trlaxlal  extension  (TXE)  tests  were 
done  on  MB  sand.  They  showed  a  slope  of  -0.88  for  the  extension  failure 
surface,  compared  with  a  slope  of  +1.39  In  trlaxlal  compression.  This 
observed  asymnetry  Is  consistent  with  test  results  on  other  sands.  The 
calculated  friction  angles  for  trlaxlal  compression  and  extension  are: 


*  2(arc  tan 


f3  -  Me 
1  *  2Me 


45* )  =  31.0° 


Trlaxlal  behavior  for  RB  sand  is  shown  In  Figure  8.  For  this 
material  the  US  and  DS  responses  indicate  slightly  different  apparent 
effective  stress  failure  surfaces,  as  seen  In  Figure  8B.  The  dry 
specimens  fall  along  a  linear  failure  surface  which  displays  a  small 
amount  of  cohesion  (0.6  MPa  Intercept  In  Figure  8B,  b'),  and  a  slope  of 
M  >  1.09.  The  calculated  friction  angle  and  cohesion  for  the  UD  tests  on 
RB  sand  are: 

„  /3  ♦  2M 

*c  «  2(arc  tan  y  3  _  c  -  45°)  =  27.5° 
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Confining  Pressure 


Figure  8.  Representative  TX  Stress-Strain  and  Effective  Stress  Path  Relations  for  RB  Sand 


c'  =  - —  -3t>l  ■  ■  =  0.29  MPa 

2V0  ♦  2MC)(3  -  Mc) 

However,  the  saturated  effective  stress  failure  envelope  is  slightly 
curved  and  lies  above  the  UD  envelope.  The  undrained  TXC  behavior  of  RB 
"sand  is  similar  to  that  of  MB  sand,  but  the  undrained  stress  paths  travel 
much  further  up  the  effective  stress  failure  envelope,  indicating  a 
greater  dilatant  tendency  (Figure  8B).  This  can  also  be  seen  in 
Figure  9.  The  relatively  large  drop  in  excess  pore  pressure  associated 
with  yielding  in  RB  sand  causes  an  increased  shear  stress  capacity.  The 
pore  pressure  Increases  upon  unloading,  which  means  the  effective  stress 
is  dropping  faster  than  total  stress  because  the  sample  is  trying  to 
compress. 

Figure  10  shows  the  trl axial  shear  stress- strain  response  for  dry  RV 
alluvium  at  several  confining  pressures,  all  considerably  higher  than 
those  used  In  the  RB  and  MB  laboratory  tests.  For  practical  purposes 
these  curves  can  be  considered  to  be  drained.  In  the  sense  that  there  was 
no  excess  pore  pressure.  The  qualitative  TX  response  of  RV  alluvium  is 
typical  of  many  dry  sands,  including  MB  and  RB  sands. 

The  observed  behavior  of  MB,  RB  and  RV  sands  under  laboratory 
conditions  can  be  summarized  as  follows: 

a)  loading  compressibility  is  nonlinear; 

b)  a  relatively  small  amount  of  strain  is  recovered  upon  unloading, 
resulting  in  substantial  net  compaction; 

c)  assumed  elastic  relationships  are  not  linear,  and  display  large 
discontinuities  at  the  transition  from  loading  to  unloading; 

d)  shear  stress- strain  response  is  gradually  softening  in  nature. 
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Axial  Strain  [%) 


Figure  10.  Typical  Stress-Strain  Response  for  Dry  RV  Alluvium  (Ref 


and  large  strains  ( >1 5  percent)  may  be  reached  before  actual 
shear  failure  occurs; 

e)  effective  stress  controls  shearing  behavior; 

f)  the  failure  surface  Is  not  symmetric  about  the  p-axls  in  the 
trl axial  plane  [o' ^  *  a'j),  hut  the  calculated  friction 
angles  In  compression  and  extension  are  approximately  equal. 
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4.0  MODELING  THE  LABORATORY  BEHAVIOR  OF  SAND 

4.1  Constitutive  Relationships 

Modeling  the  laboratory  stress-strain  response  of  any  material  is 
essentially  curve- fitting.  Curves  generated  In  the  laboratory  are  used  to 
define  the  parameters  for  a  model  of  interest.  The  model  is  then  used  to 
predict  stress-strain  response  for  conditions  not  created  in  the 
laboratory.  One  key  to  developing  a  model  which  can  accurately  predict 
stress-strain  behavior  over  a  wide  range  of  load  situations  is  using  all 
pertinent  information  from  the  available  test  results.  More  important 
than  the  number  of  model  parameters  (although  it  is  desirable  to  minimize 
that  number)  is  the  degree  to  which  all  significant  stress-strain  behavior 
mechanisms  have  been  accounted  for. 

The  primary  behaviors  studied  in  the  WES  testing  program  were: 
uniaxial  compressibility  (loading  shape  and  compaction),  shear  response 
during  uniaxial  loading  and  unloading  ( UX  stress  path),  and  triaxial  shear 
stress-strain  behavior  including  shear  failure.  Four  constitutive  models 
were  chosen  to  describe  these  behaviors  and  to  illustrate  varying  degrees 
of  thoroughness  in  utilizing  laboratory  data:  an  elastic-perfectly 
plastic  model  (ELPLA) ,  the  AFWL  Engineering  model  ( AFWL ) ,  the  Weidlinger 
effective  stress  cap  model  (CAP),  and  a  typerbolic  curve-fit  (HYPER). 

These  models  are  described  in  Appendix  B. 

4.2  Model  Data  Comparisons 

Parameters  were  selected  for  each  model  based  on  the  above 
representative  data  for  MB,  RB,  and  RV  sands.  The  model  parameters  were 
adjusted  until  they  yielded  results  which  came  within  measured  data 
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scatter,  at  least  for  those  responses  which  the  models  can  predict.  In 
some  cases,  a  model  was  developed  for  each  material  test  condition: 
undrained  dry  (UD),  drained  saturated  (DS),  and  undrained  saturated  (US). 
The  undrained  saturated  response  of  all  three  sands  is  essentially 
linearly  elastic  In  uniaxial  (confined)  compression,  with  a  von  Mises 
failure  surface  In  shear,  and  only  the  effective  stress  cap  model  was 
formulated  for  US  response.  Pore  pressure  generation  during  triaxial 
shear  was  then  predicted  with  this  model.  Specific  values  for  the  model 
parameters  are  given  In  Appendix  B. 

Figures  11  through  16  show  comparisons  between  calculated  and 
measured  responses  for  MB  sand  and  RB  sand.  RV  alluvium  results  are  vs.own 
In  Figure  17.  The  curves  labeled  "Calculation  vs.  Data"  represent  the 
integrated  differences  between  the  calculated  model  response  and  the 
measured  soil  behavior,  normalized  to  the  total  area  under  the 
representative  data  curve.  This  area  difference  is  then  divided  by  one 
hundred  to  yield  a  percent.  The  residual  has  been  calculated  for  loading 
only.  A  calculation  which  had  a  constant  area  difference  of  zero  would 
lie  exactly  along  the  data  curve.  The  rate  at  which  the  residual  deviates 
from  zero  Indicates  the  quality  of  the  model  response. 

Figure  18  shows  how  the  hyperbolic  curve-fitting  technique  may  be 
used  to  model  triaxial  stress- strain  response  at  varying  confining 
pressures.  An  example  of  results  using  the  effective  stress  cap  model  to 
simulated  undrained  behavior  Is  shown  in  Figure  19. 

4.3  Discussion  of  Comparisons 

The  elastic-plastic  model  does  fairly  well  in  defining  two  primary 
material  behaviors:  loading  compressibility  and  limited  shearing 
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Figure  12.  Drained  Saturated  MB  Sand  UX  Model -Data  Comparisons 
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Figure  16  RB  Sand  TX  Model -Data  Comparisons  (confining  pressure  =  3.45  MPa) 
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resistance.  In  the  range  of  axial  stress  used  to  test  the  sands,  a  linear 
approximation  to  the  UX  loading  curve  Is  satisfactory,  although  the 
different  sands  showed  different  degrees  of  nonlinearity.  The  materials 
tested  were  much  stlffer  on  unloading  than  on  loading,  and  recovered  only 
about  twenty  percent  of  accumulated  axial  strain.  The  elastic-plastic 
model  predicts  no  compaction  and  a  constant  for  both  loading  and  unloading 
modulus,  which  Is  a  substantial  modeling  deficiency  for  sands. 

The  AFWL  Engineering  model  predicts  nonlinear  loading  compressibility 
quite  well,  and  is  limited  only  by  the  chosen  number  of  linear  segments. 
The  model  uses  a  constant  modulus  for  unloading,  and  predicts  net 
compaction  fairly  accurately.  On  unloading  from  low  stress  levels  the 
sands  tend  to  recover  strain  at  a  greater  rate  than  from  higher  stress 
levels,  but  this  Is  not  well  modeled  by  the  AFWL  model.  This  can  be 
remedied  by  Including  a  multilinear  unload-reload  relationship,  as  shown 
In  Figure  20.  Since  Poisson's  Ratio  in  the  AFWL  model  may  be  different 
for  loading  and  unloading,  the  observed  UX  stress  path  may  be 
approximately  (l.e.,  linearly)  matched.  The  AFWL  model  triaxlal  shear 
behavior  Is  essentially  the  same  as  for  the  elastic-plastic  model.  The 
ultimate  shear  strength  Is  well  predicted  by  the  linear  failure  surface, 
but  the  shape  of  the  TX  stress-strain  curve  Is  not.  The  shear  modulus  is 
directly  proportional  to  the  bulk  modulus,  and  so  the  model  may  even 
predict  a  stiffening  shear  modulus  If  a  Kydrostat  breakpoint  is 
encountered  during  triaxlal  shear.  The  poor  prediction  of  shear 
stress-strain  behavior  Is  a  weakness  of  the  AFWL  model. 

The  cap  model  (as  Installed  In  the  SEM)  predicts  only  a  concave  (to 
the  stress  axis)  UX  loading  compressibility  curve.  Thus,  changes  In 
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Figure  20.  AFWL  Engineering  Model  with  Multilinear 
Reload  Behavior 
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curvature  In  the  UX  curve  are  not  modeled.  Most  of  the  sand  data  is 
consistently  convex,  so  the  cap  model  does  as  good  a  job  as  does  the  AFWL 
model.  The  nonlinear  nature  of  the  UX  unloading  behavior  is  handled  much 
better  by  the  cap  model,  as  seen  in  the  model-data  comparisons  (e.g. 

-Fig.  13).  The  UX  stress  path  is  somewhat  hard  to  control,  but  can  be  made 
to  follow  the  trend  of  the  data.  Since  the  cap  model  parameters  are  fit 
to  follow  TX  stress-strain  curvature,  the  model  is  considerably  more 
accurate  In  predicting  these  results  than  are  either  the  elastic-plastic 
or  AFWL  models.  The  cap  model  cannot  predict  a  post  peak  softening 
behavior,  as  was  observed  in  drained  saturated  RB  sand  (Fig.  16). 

The  cap  model  was  used  to  predict  undrained  behavior,  since  it  can 
calculate  both  total  and  effective  stress.  The  primary  cap  modeling 
deficiency  In  this  exercise  (as  seen  In  Figure  19  for  RB  sand)  was  its 
inability  to  follow  the  effective  stress  failure  envelope  upwards  after 
Initially  reaching  It.  In  the  actual  material  response,  a  significant 
post  "failure"  Increase  in  shear  strength  was  observed,  due  to  the 
sample's  tendency  to  dilate.  The  model  does  not  consider  this  behavior, 
and  therefore  did  not  do  well  In  predicting  stress-strain  relations,  shear 
strength,  or  pore  pressure  response.  A  total  stress  approach  with  a  von 
Mlses  type  failure  surface  fit  to  ultimate  strength  would  be  more 
accurate,  but  fundamentally  lacking,  since  effective  stress  actually 
controls  shear  behavior.  The  problem  lies  in  predicting  excess  pore 
pressure,  l.e.,  In  predicting  the  coupling  between  shear  strain  and  the 
tendency  of  the  soil  skeleton  to  change  volume. 
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5.0  DYNAMIC  INSITU  BEHAVIOR  OF  SAND 

5.1  Insltu  Test  Description 

The  insltu  behavior  of  sand  will  be  discussed  here  with  reference  to 
’two  types  of  Insltu  dynamic  material  property  tests:  the  Cylindrical 
Insltu  Test  (CIST)  and  the  Dynamic  Insitu  Compressibility  Test  (DISC 
Test).  These  tests  have  each  been  performed  at  the  Ralston  Valley 
Geotechnical  Study  Site  In  Nevada.  The  laboratory  properties  of  RV 
alluvium  are  available  (Ref.  11)  and  have  been  discussed  previously. 

CIST  22  (Ref.  12)  fielded  gages  on  several  redundant  radials  at  a 
single  depth  of  4.72  m.  Accelerometers  measured  radial  and  vertical 
motions,  and  flat  pack  stress  gages  measured  stress  time  histories.  DISC 
Tests  I  and  II  (Refs.  13  and  14)  were  performed  at  approximately  the  same 
site  as  CIST  22.  Their  purpose  was  to  define  insltu  vertical  uniaxial 
strain  compressibility  to  a  depth  of  about  fifteen  meters.  Approximately 
thirty  ground  motion  gages  were  placed  for  each  test  (vertical  and 
horizontal)  along  with  10  alrblast  gages.  (Before  edge  effects  appear, 
the  tests  maintain  a  condition  of  essentially  uniaxial  strain.) 

Figure  21  shows  the  experimental  configuration  for  CIST  22,  and  Figure  22 
shows  the  experimental  configuration  for  DISC  Test  II. 

5.2  Insltu  Test  Results 

Several  Important  features  of  CIST  and  DISC  Test  data  help  define 
Insitu  material  properties: 

1)  Attenuation  of  peak  motions  (acceleration,  velocity,  and 
displacement; 

11)  Arrival  times  of  first  motion  and  peaks; 
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Ill)  Waveform  shape  (rise  time,  post  peak  rate  of  decay,  and  number 
of  cycles  of  motion); 

1v)  Two-dimensional  effects  as  reflected  In  the  waveforms. 

Figure  23  shows  composite  velocity  waveforms  and  Figure  24  shows  how 
j>eak  velocity  and  arrival  times  varied  with  range  in  CIST  22.  The  data  Is 
consistent  with  other  tests  In  dry  sand,  although  the  quality  (based  on 
redundancy)  Is  not  good.  Definition  of  the  actual  cavity  pressure  was  not 
totally  successful,  but  some  gages  did  survive  and  a  reasonable  estimate 
can  be  made  based  on  the  gages  which  did  survive. 

The  primary  (but  not  the  only)  difference  between  a  DISC  Test  and  a 
CIST  Is  the  direction  of  motion.  A  CIST  is  designed  to  test  each 
individual  soil  1*yer  by  Inducing  radial  displacement.  A  DISC  Test 
Includes  the  effects  of  reflection  and  refraction  due  to  layering  by 
Inducing  vertical  displacement.  Waveforms  are  shown  for  DISC  Test  II  in 
Figure  25.  Note  the  effect  of  layering  (l.e.  reflections)  and  explosive 
cavity  venting  through  the  berm,  which  give  the  waveforms  more  character 
than  shown  by  CIST  waveforms.  Attenuation  and  arrival  times  for  DISC 
Test  II  are  shown  In  Figure  26. 

Both  the  CIST  and  DISC  Test  results  are  Influenced  by  two-dimensional 
effects.  Relief  waves  propagate  Into  the  fields  of  recorded  motions  from 
the  edge  of  the  explosive  cavity,  and  the  free  surface  affects  propagating 
waves  In  both  test  geometries. 
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B.  Arrival  Times 

Figure  24.  CIST  22  Observed  Attenuation  and  Arrival  Times 
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Figure  26.  DISC  Test  II  Observed  Attenuation  and  Arrival  Times 
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6.0  MODELING  THE  INSITU  BEHAVIOR  OF  SAND 


6.1  Finite  Difference  Calculations 

The  CIST  and  DISC  Test  events  In  Ralston  Valley  alluvium  were 
simulated  with  a  series  of  finite  difference  calculations.  The  majority 
of  calculations  were  either  one- dimensional ,  axisymmetrlc  or  plane- 
strain.  The  one-dimensional  code  Installed  In  the  Soil  Element  Model  (see 
Section  2.1)  was  used.  In  addition,  a  two-dimensional  plane-strain 
calculation  was  done  for  DISC  Test  II  using  STEALTH-2D  (Ref.  15).  This 
illustrated  the  two-dimensional  nature  of  motions  generated  in  the  DISC 
Test.  The  boundary  conditions  and  grid  definition  for  the  one-dimensional 
calculations  are  shown  In  Figure  27.  Edge  and  bottom  boundary  conditions 
were  not  Important,  because  the  grid  was  large  enough  to  provide  ample 
simulation  time  free  of  reflections. 

Several  constitutive  relationships  were  used  to  calculate  each 
event:  a  11  near- el as tic  relationship,  a  simple  elastic-plastic  model,  the 
AFWL  Engineering  model,  and  the  cap  model.  The  model  parameters  were 
those  defined  by  laboratory  test  results  (see  values  given  In 
Appendix  B).  Additional  calculations  were  done  with  revised  model 
parameters  to  better  match  the  measured  waveforms.  Table  3  provides  a 
summary  of  the  calculations  performed. 

6.2  Calculatlonal  Results 

Results  generated  for  each  1-D  calculation  consist  of: 

a.  Time-histories  of  stress,  velocity,  displacement,  and  any  number 
of  additional  arbitrary  quantities  (e.g.  cap  position); 

b.  Snapshots  of  stresses  and  strains  In  the  grid  at  various  times; 
and 
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B.  DISC  Test  II  Calculation 


Figure  27. 


One-Dimensional  Calculational  Setups  for 
CIST  22  and  DISC  Test  II 


Table  3.  Finite-Difference  Calculation  Matrix 


tvent 

Calc  No. 

Material  Model 

Notes 

CIST  22 

C22-1 

Linear  Elastic 

No  tensile  failure 

C22-2 

Elastic-Plastic 

C22-3 

AFWL 

Laboratory  based  parameters 

C22-4 

CAP 

Laboratory  based  parameters 

C22-5 

AFWL 

Insltu  based  parameters 

DISC  Test  II 

DT-1 

Linear  Elastic 

DT-2 

Elastic-Plastic 

DT-3 

AFWL 

Laboratory  based  parameters 

DT-4 

CAP 

High-pressure  laboratory  fit 

DT-5 

CAP 

Low-pressure  laboratory  fit 

DT-6 

AFWL 

2-D,  laboratory  parameters 

c.  SEK  constitutive  model  plots  at  one  target  point. 

Appendix  C  contains  an  example  a  set  of  complete  results  for  a 
calculation.  Three  aspects  of  the  calculations  will  be  compared  here: 
waveform  shape,  attenuation  of  peak  values,  and  arrival  times. 

The  elastic  calculation  of  CIST  22  did  not  produce  waveforms  which 
resembled  the  observed  motion,  as  shown  In  Figure  28.  The  remainder  of 
the  models,  however,  did  match  the  general  character  of  the  data,  as  seen 
In  Figure  29.  Attenuation  comparisons  between  the  various  models  and  the 
CIST  data  are  shown  In  Figure  31A. 

The  character  of  motion  produced  In  the  DISC  Test  calculations  was 
quite  similar  for  all  models.  Figure  30  compares  calculated  waveforms 
with  data  near  the  center  of  the  explosive  bed.  (The  elastic  model  Is  not 
Included,  as  It  overlaid  the  elastic-plastic  calculation.)  Attenuation  is 
shown  In  Figure  31B.  Time  of  arrival  of  the  first  signal  and  velocity 
peaks  are  shown  for  CIST  and  DISC  Test  II  in  Figure  32. 

Results  from  the  two-dimensional  DISC  Test  calculation  are  summarized 
In  Figure  33.  The  figure  shows  how  vertical  and  horizontal  velocity 
waveforms  vary  throughout  the  test  bed. 

Some  general  conclusions  can  be  drawn  from  the  above  calculations: 
a)  Necessary  model  complexity  varies  with  calculatlonal  geometry. 

The  elastic,  elastic-plastic,  AFWL,  and  cap  models  all  produced  similar 
results  for  the  DISC  Test.  This  Is  because  motion  is  uniaxial,  and  the 
most  Important  material  response  modeled  Is  uniaxial  compression,  which 
all  four  models  do  fairly  well  within  their  bounds.  Attenuation  Is  a  main 
difference  between  models.  Note  that  an  elastic  model  should 
theoretically  not  attenuate  motion  or  stress  at  all.  However,  artificial 
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Figure  28.  CIST  22  Elastic  Calculation  Compared  with  Data 
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Figure  31.  Comparison  of  Calculated  and  Measured  Attenuation  Rates 
for  CIST  22  and  DISC  Test  II 
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Figure  32.  Comparison  of  Calculated  and  Measured  Times  of  Arrival 
for  CIST  22  and  DISC  Test  II 


Range  (m) 


Figure  33.  DISC  Test  II  Calculated  (2D  STEALTH)  Velocity  Waveforms 


viscosity  reduces  peaks,  and  this  effect  was  seen  In  Figures  316.  In  the 
CIST  calculations  the  results  change  markedly  In  going  from  an  elastic 
model  to  an  elastic-plastic  model  with  shear  and  tensile  failure.  This  Is 
due  to  the  cylindrical  symmetry  of  the  problem.  Stress  relief  In  the  hoop 
direction  generates  a  more  complex  stress/straln  field  than  In  the 
uniaxial  compression  case,  even  In  one  dimensional  (radial)  wave 
propagation.  In  order  of  Importance,  the  critical  features  of  material 
response  In  the  CIST  calculations  appear  to  be: 

1)  Volumetric  (compresslonal)  behavior; 

11)  Shear  failure  and  post  failure  behavior;  and 

111)  Tensile  behavior. 

Cylindrical  calculations  therefore  tend  to  reveal  much  more  than 
planar  calculations  about  the  constitutive  model  employed. 

b)  Model  parameters  should  be  based  on  laboratory  data  at  the  same 
level  of  stress  expected  Insltu.  The  cap  model  (which  predicts  Insltu 
motions  considerably  larger  than  observed)  was  fit  at  the  70  lf>a  level  in 
the  uniaxial  test,  while  actual  stress  levels  were  on  the  order  of 

5-10  H>a.  At  lower  stresses  the  model  Is  clearly  too  soft,  but  can  easily 
be  adjusted  to  better  fit  low  stress  level  laboratory  data.  If  a  wide 
range  of  peak  stress  is  expected  In  a  particular  problem,  the  model  should 
match  the  data  at  all  stress  levels  below  the  maximum.  Although  the 
models  could  be  modified  to  facilitate  curve  fitting,  detailed  matching 
over  a  wide  stress  range  was  not  possible  with  the  Implemented  versions  of 
the  AFWL  and  cap  models. 

c)  Two-dimensional  effects  can  Influence  Insltu  model  estimation. 
Based  on  waveform  matching  In  one-dimensional  calculations,  an  insitu- 
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based  model  can  be  selected.  Inherent  In  the  measured  waveforms  are 
relief  and  other  two-dimensional  effects,  which  must  be  recognized  during 
modeling.  Otherwise,  the  model  may  be  biased  by  unloading  effects,  and 
may  therefore  significantly  underpredict  motion. 
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7.0  DISCUSSION  OF  SAND  MODELING 

7.1  Additional  Laboratory  Behavior  of  Sand 

More  Is  known  about  the  laboratory  behavior  of  cohesionless  soil  than 
Is  reflected  In  the  data  available  for  this  study.  The  following  is  a 
summary  of  points  which  would  be  important  in  developing  a  complete 
constitutive  relationship. 

7.1.1  Multiphase  Effects.  Soil  is  a  two-or  three-phase  medium,with 
a  solid  soil  skeleton  and  water,  air,  or  both  filling  the  voids.  The  rate 
of  pore  fluid  migration  through  the  skeleton  has  a  significant  effect  on 
effective  stress,  and  therefore  on  behavior  under  both  static  and  dynamic 
loads.  The  permeability  of  sand  depends  to  a  large  degree  on  the  amount 
of  fines  present,  and  Is  usually  quite  high  compared  with  that  for  silts 
and  clays.  Depending  on  the  type  and  rate  of  loading,  pore  fluid 
migration  In  sand  may  or  may  not  occur.  Lacking  data.  It  Is  necessary  to 
estimate  the  permeability  of  a  sand  based  on  void  ratio  and  percent  of 
fines. 

When  a  soil  Is  partially  saturated,  its  behavior  Is  considerably 
different  from  that  In  either  the  dry  or  the  fully  saturated  condition. 
Pore  pressure  parameters  may  be  defined  to  describe  partially  saturated 
behavior  as  done  by  Skempton  In  Reference  16.  Under  dynamic  loads,  pore 
air  may  enter  Into  solution,  and  the  soil  may  become  fully  saturated.  The 
details  of  this  phenomenon  need  to  be  quantified  for  modeling  purposes. 

7.1.2  Anlstropy  and  Work  Hardening.  Soil  in  Its  natural  state  may 
be  Inherently  anisotropic,  or  It  may  become  anisotropic  due  to  application 
of  stress.  A  grain  matrix  with  particles  oriented  In  a  preferred 
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direction  will  complicate  definition  of  a  constitutive  relationship 
because  complete  symmetry  can  no  longer  be  assumed. 

As  a  soil  deforms.  Its  strength  characteristics  often  change.  Two 
types  of  hardening/softening  behavior  have  been  postulated  for  soil: 
dsotroplc  and  kinematic.  Figure  34  illustrates  the  basic  concept  behind 
each  of  these  phenomena.  For  an  isotropically  hardening  material,  the 
loading  surface  expands  uniformly  about  the  origin  in  stress  space  and 
maintains  the  same  shape,  center,  and  orientation  as  the  yield  surface. 

In  kinematic  hardening,  plastic  deformation  causes  the  loading  surface  to 
translate  as  a  rigid  body  In  stress  space.  Changes  in  properties  as  a 
result  of  loading  are  Important  for  long  duration  and  multi-directional 
loads,  and  for  cases  of  multiple  loading. 

7.1.3  Strains.  It  has  been  well  established  that  sands  often  change 
(or  tend  to  change)  volume  as  a  result  of  shear  deformation.  This  Is 
known  as  shear-volume  coupling,  or  dllatancy.  This  coupling  Is  especially 
Important  In  determining  pore  pressure  response  In  saturated  soils.  Some 
evidence  of  this  was  noted  in  the  tri axial  behavior  of  MB  and  RB  sand,  but 
the  exact  relationship  for  coupled  shear-volume  behavior  was  not  well 
defined. 

Another  aspect  of  strain  behavior  not  explicitly  covered  thus  far  is 
the  definition  of  plastic  yield  and  potential  functions.  It  has  been 
postulated  that  soil  deforms  In  the  direction  of  the  stress  Increment 
under  small  stress  Increments,  and  In  the  direction  of  total  stress  (and 
not  stress  Increment)  under  large  stress  increments.  It  is  therefore 
necessary  to  define  a  plastic  potential,  in  addition  to  a  failure  or  yield 
surface,  which  defines  the  plastic  strain  behavior  of  the  material.  The 
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nature  of  the  plastic  potential  with  respect  to  the  yield  surface 
determines  whether  the  flow  rule  is  associative  or  nonassoclative. 

7.1.4  Loading  Rate  Effects.  Test  results  on  MB  and  RB  sands 
Indicate  that  properties  were  not  dependent  on  loading  rate,  for  the 
^loading  rates  employed.  It  has  been  shown  by  Jackson,  et  al.  (Ref.  17) 
that  only  submillisecond  loading  rates  influence  the  compressibility  of 
cohesionless  sands.  It  Is  assumed  that  sandy  soils  with  an  increasingly 
higher  percentage  of  fines  are  more  sensitive  to  loading  rate.  This 
relationship  needs  to  be  quantified  for  poorly  sorted  alluvium,  or  silty 
soils  such  as  playas. 

7.1.5  Cyclic  Loading  Effects.  The  laboratory  tests  and  calculations 
In  this  study  dealt  with  explosive  or  Impulsive  loadings,  which  generate 
motions  having  only  a  few  cycles  at  most.  However,  a  substantial  amount 
of  research  has  been  done  to  Investigate  cyclic  loading  effects  on  soils, 
because  of  their  Importance  In  other  dynamic  loading  phenomena.  Some  of 
the  more  Important  aspects  of  cyclic  behavior  are: 

a)  Energy  Dissipation.  Each  cycle  of  loading-unloading  with 
hysteresis  (even  at  low  stress  levels)  represents  a  loss  of 
energy  through  material  damping.  Damping  characteristics  are 
Important  for  problems  such  as  foundation  vibration. 

b)  Liquefaction.  The  tendency  for  changes  In  pore  volume  to  occur 
during  cyclic  loading  will  cause  excess  pore  water  pressure 
buildup.  This  results  In  a  drop  In  effective  stress,  which  may 
cause  shear  failure. 

c)  Modulus  Degradation.  As  strains  accumulate  under  many  cycles  of 
load,  the  shear  modulus  may  decrease,  resulting  In  an  observed 
softening  behavior. 


7.1.6  High  Pressure/Temperature  Behavior.  At  very  high  pressures 
and  temperatures  the  solid  constituents  of  soil  may  undergo  phase 
changes.  These  changes  have  Important  energy  balance  consequences  for  the 
soil  system.  The  nature  of  these  phenomena  is  typically  explored  through 
/Iyer  plate  and  gas  gun  experiments. 

7.2  Additional  Aspects  of  the  Insitu  Behavior  of  Sand 

The  calculations  presented  in  this  report  address  a  specific  aspect 
of  dynamic  soil  behavior,  insitu  material  property  definition  through  high 
explosive  tests.  There  are  mary  other  dynamic  loading  situations  which 
involve  additional  responses  fundamental  to  soils.  They  can  be  grouped 
into  the  following  categories: 

a)  Explosive  Loading.  These  are  high  amplitude  impulsive  loads  with 
relatively  few  cycles  of  motion.  Several  types  of  explosive  load 
configurations  are  Illustrated  in  Figure  35.  Point  sources  may 
be  either  nuclear  or  high  explosive,  and  vary  in  position  with 
respect  to  the  ground  surface.  HEST  and  DIHEST  type  loads  are 
used  to  simulate  nuclear  environments,  or  to  determine  Insitu 
material  properties.  There  is  a  large  amount  of  data  on 
explosive  effects  In  sandy  soils,  as  shown  In  Table  4.  The  bulk 
of  the  data  base  comes  from  a  few  test  sites,  and  is  mostly  for 
dry  alluvium. 

b)  Impact.  Impact  loads  are  similar  to  explosive  loads,  in  that 
they  are  very  Impulsive  and  produce  few  cycles  of  motion. 

c)  Earthquake.  These  loads  usually  have  many  cycles,  and  durations 
on  the  order  of  seconds.  Transmission  Is  typically  in  bedrock, 
and  energy  propagates  to  the  surface  through  the  soil.  A  common 


66 


D.*v  o,Itj 


un-ace 


o 


Planar  Array 

Cylindrical  -  Row  (DIHEST) 
Cylindrical  -  Single  (CIST) 


C.  Direct  Induced 


Figure  35.  Types  of  Explosive  Loading 
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assumption  Is  that  the  waves  In  the  soil  are  vertically 
propagating  shear  waves. 

d)  Wave-Induced.  Water  waves  produce  oscillatory  loads  In 
underlying  sediments. 

e)  Vehicular.  The  wheels  or  tracks  of  a  vehicle  moving  across  the 
ground  surface  produce  transient  loads  of  a  somewhat  complicated 
nature. 

f)  Foundation  Vibration.  These  loads  are  transmitted  from  a  moving 
foundation  Into  the  surrounding  soil.  Excitation  may  be  due  to 
machine  vibrations,  wind  loads,  impact,  etc.  Thus  these  loads 
can  be  either  transient  or  steady  state  In  nature,  and  can  vary 
widely  In  magnitude  and  frequency. 

Each  type  of  dynamic  loading  may  produce  a  different  type  of 
response.  Several  Important  types  of  dynamic  soil  behavior  include: 

a)  Transient  Waves.  Energy  transmission  through  soil  can  be 
characterized  by  wave  and  phase  velocity,  peak  stress,  peak 
motion  (acceleration,  velocity,  displacement),  attenuation  of 
peaks  (energy  dissipation),  frequency  content,  rise  time  and 
duration. 

b)  Steady-State  Waves.  A  stationary  wave  pattern  requires  a  steady- 
state  source,  due  to  the  dissipative  nature  of  soil. 

c)  Cratering.  The  separation  and  ejection  of  material  around  a 
local  source  of  energy. 

d)  Spall.  Soil  may  become  distended  (loss  of  particle-particle 
contact);  also  called  lofting. 

e)  Creep.  Long  term,  time- dependent  response. 
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f)  Liquefaction.  Loss  of  strength  due  to  sudden  drop  In  effective 
stress. 

Table  5  shows  the  phenomena  which  are  particularly  Important  for  each 
category  of  dynamic  loading. 

The  dynamic  phenomena  mentioned  above  are  manifestations  of  various 

T 

aspects  of  soil  behavior.  Thus,  on  a  fundamental  level,  soil  behavior  can 
be  tied  to  specific  observable  phenomena,  and  this  is  done  In  Table  6. 

This  table  allows  a  material  modeler  to  choose  the  behaviors  which  must  be 
accounted  for  when  a  particular  response  Is  being  calculated. 


a 

Table  5.  Important  Response  Phenomena  for  Dynamic  Loadings 


Important  Response 
Phenomena 

Transient  Waves 

Steady-State  Waves 

Crateri ng 

Spall 

Creep 

Liquefaction 


Foundation  Vibration 


Table  6.  Controlling  Soil  Behavior  for  Various 
Dynamic  Response  Phenomena  in  Sands 


T 

Controlling  Soil  Behavior 

— 

Dynamic  Response  Phenomena 

Transient  Waves 

i - 

Steady  State  Waves 
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X 
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X 
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X 

Pressure  Dependent  Shear  Strength 
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II 
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n 

X 
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II 

X 

X 

X 

X 

Pore  Fluid  Movement 

9 

X 

Minimal  Tensile  Strength 

II 

X 

X 

Phase  Changes 

X 

8.0  TESTING  CONSIDERATIONS 

8.1  Ideal  Testing  Program 

The  purpose  of  a  soil  stress-strain  test  is  to  obtain  data  with  which 
to  predict  the  load-deformation  response  of  a  soil  mass,  for  a  particular 
field  loading  condition.  In  the  context  of  the  present  research,  the 
field  loading  of  particular  interest  is  an  explosion,  either  nuclear  or 
chemical,  detonated  above,  at,  or  below  the  ground  surface.  The  problem 
is  how  to  make  the  prediction,  i.e. ,  what  stress-strain  data  the  soil  test 
should  yield,  and  how  that  stress-strain  data  should  be  used  to  calculate 
load-deformation  response. 

There  are  two  principal  philosophies  or  methods  of  soil  stress-strain 
testing.  The  first  is  the  Load  Path  philosophy.  Under  this  philosophy 
the  soil  stress-strain  test  is  designed  to  duplicate  a  representative 
field  loading  (stress  or  strain)  path.  The  test  data  is  then  used 
*i;rertly  in  the  field  load- deformation  prediction.  An  early  example  of 
this  approach  is  use  of  confined  compression  (oedometer)  data  to  predict 
foundation  settlement.  The  two  main  problems  with  this  method  are: 
knowing  what  loading  paths  to  use  In  the  test;  and  for  all  but  the 
simplest  loading  paths,  constructing  a  test  device  capable  of  duplicating 
the  desired  loading  paths. 

The  explosive  loading  problem  Involves  a  large  variation  in  loading 
amplitudes,  the  Interaction  of  propagating  waves  with  materials  of 
different  properties  and  the  free  surface,  the  creation  of  new  waves  as  a 
result  of  these  Interactions  and  the  changes  In  the  material  states  as 
these  waves  Interact.  Simplifying  these  complex  conditions  to  a  few 


states  of  stress  or  strain  has  proven  very  difficult.  Trullis  (Ref.  9) 
has  suggested  that  for  a  burled  burst  In  a  homogeneous  media  the  strain 
paths  are  quite  similar  and  Independent  of  the  properties  of  the  media. 
Experiments  to  demonstrate  the  accuracy  of  this  statement  based  on 
theoretical  calculations  have,  to  date,  been  unsuccessful.  The  strain 

•-r 

paths  calculated  for  surface  burst  are  more  complex  and  vary  more, 
however,  Trullls  still  maintains  that  they  may  be  generalized. 

Experimental  evaluation  for  this  case  Is  beyond  the  state-of-the-art  In 
ground  motion  instrumentations.  This  analysis  further  requires  that  the 
soil  sample  be  subjected  to  a  tensile  strain  in  a  principal  direction, 
which  presents  serious  problems  In  developing  a  laboratory  device  to 
subject  a  sample  to  the  required  strain  paths  and  allow  measurement  of  the 
associated  stresses. 

The  second  philosophy  Is  the  Constitutive  Model  philosophy.  Under 
this  philosophy  one  assumes  the  existence  of  a  mathematical  constitutive 
model,  sufficiently  general  and  accurate  to  correctly  predict  soil  stress- 
strain  response  for  loading  paths  other  than  those  used  to  formulate  the 
model  and  evaluate  Its  parameters.  The  constitutive  model  need  not  be  the 
"one  and  only  true  model",  but  it  must  yield  accurate  predictions  of  load- 
deformation  response  for  field  loading  cases  of  interest.  The  main 
elements  In  this  method  are:  establishing  what  stress-strain  mechanisms 
to  model  (and  which  ones  to  Ignore),  selecting  a  satisfactory  model, 
deciding  what  tests  are  needed  to  evaluate  the  model  parameters, 
constructing  a  test  device  capable  of  performing  the  required  tests,  and 
finally  solving  the  field  loading  boundary  value/propagation  problem  using 
the  selected  constitutive  model. 
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The  Ideal  soil  stress-strain  testing  philosophy  is  a  compromise 
between  (or  a  combination  of)  the  above  two  philosophies.  It  involves 
formulating  a  constitutive  model  based  on  test  load  paths  as  close  as 
possible  to  the  field  loading  paths.  This  approach  combines  the  realism 
»of  the  Load  Path  Method  with  the  generality  of  the  Constitutive  Model 
Method. 

The  main  decisions  to  be  made  in  designing  the  above  ideal  laboratory 
soil  stress-strain  testing  program  for  a  given  field  (in  this  case 
explosive)  loading  situation  involve: 

a)  estimating  the  field  loading  path(s)  and  rates; 

b)  selecting  the  constitutive  model; 

c)  defining  the  test  loading  path(s)  and  rates; 

d)  stress  versus  strain  control; 

e)  effective  versus  total  stress; 

f)  sample  shape; 

g)  sample  size; 

h)  number  of  tests;  and 

i)  treatment  of  data  scatter. 

For  a  typical  field  explosive  loading  situation,  the  soil  is  initially  in 
what  approximates  a  KQ  condition,  with  an  OCR  _>  1.  The  airslap- induced 
motion  involves  primarily  vertical  compression;  and  the  upstream-induced 
motion  Involves  primarily  radial  compression,  outward  radial  displacement, 
and  vertical  compression.  The  outward  radial  displacement  tends  to 
produce  circumferential  tension;  at  least  it  relieves  some  of  the  Poisson 
compression  caused  by  the  radial  and  vertical  compression.  In  addition, 
the  combination  of  airslap  and  upstream-induced  motion  causes  s  rotation 
of  principal  planes  about  a  circumferential  axis. 
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One  conceptually  attractive  soil  stress-strain  test  device  for  the 
above  situation  Is  a  strain-controlled  hollow  cylinder  device.  In  which 
the  Inner  and  outer  drums  are  capable  of  applying  vertical  shear  stresses 
to  the  sample's  Inside  and  outside  curved  surfaces.  A  less  Ideal  device, 
but  one  which  already  exists.  Is  a  stress-controlled  hollow  cylinder 
device,  equipped  to  apply  a  torque  along  the  cylinder  axis.  Both  hollow 
cylinder  devices  have  the  advantage  of  being  able  to  induce 
(circumferential )  tension  using  (radial)  compressive  loading,  but  also 
have  the  disadvantage  of  creating  a  nonhomogeneous  state  of  stress  and 
strain  In  the  sample.  In  addition,  the  axis  of  principal  stress  rotation 
in  the  stress-controlled  hollow  cylinder  test  Is  radial,  rather  than 
circumferential.  Both  devices  will  yield  stress-strain  and  strength  data 
for  which  the  three  principal  stresses  are  distinct.  A  rapid  method  for 
plotting  such  data  In  the  octahedral  plane  has  been  devised  by  Merkle 
(Ref.  18). 

8.2  Modeling  Requirements 

The  two  principal  requirements  of  a  soil  constitutive  model  for  field 
explosive  loading  applications  are  that  It  Include  all  Important  physical 
mechanisms,  and  be  computationally  efficient. 

The  principal  requirements  for  the  above  model  are: 

a)  behavior  Is  controlled  by  effective  stress; 

b)  considers  pore  fluid  compressibility  and  flow  characteristics; 

c)  dynamic  hydrostat  may  be  concave  to  the  strain  axis  at  low 
stress,  but  Is  convex  to  the  strain  axis  at  high  stress; 

d)  shear  (deviator)  stress-strain  curve  concave  to  the  strain  axis 
up  to  the  point  of  peak  strength,  but  may  exhibit  strain 
softening  beyond  that  point  to  approach  a  lower  residual  strength 


e)  peak  and  residual  shear  strengths  Increase  with  Increasing 
confining  pressure,  with  nonlinear  effects  often  Important  at 
very  high  confining  pressure; 

f)  very  low  tensile  strength,  with  post  tensile  failure  stress 
redistribution; 

T 

g)  Inelastic  volumetric  and  deviatoric  strains  accumulate  from  the 
onset  of  loading; 

h)  exhibits  both  volumetric  and  deviatoric  hysteresis; 

1)  exhibits  shear  strain/volume  strain  coupling  (dllatancy), 
controlled  by  previous  strain  history; 

j)  plastic  flow  rule  may  be  nonassoclatlve  (l.e.,  plastic  strain 
Increments  may  not  be  normal  to  the  yield  surface); 

k)  accomodates  a  geologically  realistic  Initial  (KQ  consolidation) 
condition,  with  an  OCR  >  1. 

None  of  the  soil  constitutive  models  yet  exercised  by  the  SEM  satisfied 
all  the  above  requirements.  However,  the  AFWL  Engineering  and  cap  models 
both  have  several  desirable  features.  Including  familiarity  within  the 
groundshock  community.  They  are  therefore  prime  candidates  for  possible 
improvement  In  the  next  phase  of  this  research,  the  objective  of  which  is 
to  develop  a  model  satisfying  as  many  of  tfe  above  requirements  as 
possible. 
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9.0  SUMMARY 

The  FYb2  phase  of  the  study  Fundamental  Properties  of  Soils  for 
Complex  Dynamic  Loadings  focused  on  dynamic  constitutive  modeling  of  sandy 
„so11s.  Four  constitutive  models  of  explosively  loaded  sand  were  studied: 

a)  elastic-plastic; 

b)  AFWL  Engineering; 

c)  cap; 

d)  hyperbolic. 

A  one-dimensional  wave  propagation  feature  was  added  to  the  Soil  Element 
Model  (SEM)  to  study  the  Influence  of  a  constitutive  model  on  wave 
propagation  effects.  Insltu  material  property  tests  CIST  22  and 
DISC  Test  II  were  modeled  In  one  and  two  dimensions,  using  the  above  four 
constitutive  models,  and  the  calculated  and  measured  results  were 
compared.  The  SEM  was  also  used  to  assess  the  performance  of  the  four 
constitutive  models  against  the  results  of  laboratory  tests  on  three  sands 
conducted  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station. 

A  principal  conclusion  to  be  drawn  from  the  above  comparisons  of 
calculated  and  measured  results  is  that  dllatancy,  the  coupling  between 
shear  and  volumetric  strain,  needs  additional  attention  In  constitutive 
modeling  of  cohesionless  soils. 

The  conceptually  Ideal  laboratory  soil  testing  device  for  explosive 
loading  applications  Is  a  strain-controlled  hollow  cylinder  device. 
Theoretically  It  has  the  capability  to  Induce  radial  compression,  hoop 
expansion  due  to  outward  radial  displacement,  vertical  compression,  and 
rotation  of  principal  planes  about  the  circumferential  axis. 
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APPENDIX  A 


WES  Laboratory  Tests 
(Refs.  5  and  6) 

The  isotropic  compression  (IC)  test  subjects  a  cylindrically  shaped 
specimen  to  an  equal  all-around  confining  pressure  while  measurements  of 
the  specimen's  height  and  diameter  changes  are  made.  The  data  are 
normally  plotted  as  pressure  versus  volumetric  strain,  the  slope  of  which 
Is  the  bulk  modulus,  K. 

The  trl axial  compression  (TXC)  test  is  conducted  after  a  desired 
confining  pressure  Is  applied  during  the  IC  test.  While  the  confining 
pressure  Is  held  constant,  axial  load  Is  Increased  and  measurements  of  the 
specimen's  height  and  diameter  changes  are  made.  The  data  can  be  plotted 
as  principal  stress  difference  versus  axial  strain,  the  slope  of  which  is 
Young's  modulus  E,  or  as  principal  stress  difference  versus  principal 
strain  difference,  the  slope  of  which  Is  twice  the  shear  modulus  G.  The 
maximum  principal  stress  difference  the  specimen  can  support  or  the 
principal  stress  difference  at  15  percent  axial  strain  during  shear 
loading  (whichever  occurs  first)  Is  defined  as  failure  and  describes  one 
point  on  a  failure  surface  under  positive  principal  stress  difference 
states  of  stress.  The  failure  surface  is  depicted  as  a  plot  of  principal 
stress  difference  versus  mean  normal  stress. 

The  triaxlal  extension  (TXE)  test  Is  also  conducted  after  a  desired 
confining  pressure  Is  applied  during  the  IC  test.  While  lateral  pressure 
Is  held  constant,  vertical  pressure  Is  decreased  and  measurements  of  the 
specimen's  height  and  diameter  changes  are  made.  As  with  the  TXC  test, 
the  data  are  plotted  as  principal  stress  difference  versus  axial  strain  or 
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as  principal  stress  difference  versus  principal  strain  difference.  The 
maximum  negative  principal  stress  difference  or  the  point  at  which  the 
material  separates  (whichever  occurs  first)  is  defined  as  failure  and 
describes  one  point  on  a  failure  surface  under  a  negative  principal  stress 
difference  state  of  stress. 

Three  types  of  uniaxial  strain  ( UX )  tests  were  conducted: 

1)  The  first  (designated  UX)  Is  conducted  by  applying  an  axial 
(vertical)  pressure  to  a  wafer-shaped  specimen  that  Is  physically 
constrained  from  deflecting  radially.  Measurements  are  made  of 
the  applied  axial  stress  and  the  specimen's  height  change.  The 
data  are  plotted  as  axial  (vertical)  stress  versus  axial 
(vertical)  strain,  the  slope  of  which  is  the  constrained  modulus 
M. 

2)  The  second  type  of  UX  test  (designated  UX/KQ)  is  conducted  by 
applying  radial  pressure  to  a  cyllndrically  shaped  specimen  until 
a  slight  Inward  movement  of  the  diameter  is  detected.  Axial  load 
Is  then  applied  until  the  specimen  returns  to  its  original  radial 
position  (zero  radial  strain).  This  process  Is  repeated 
throughout  the  loading  and  unloading.  As  in  the  UX  test,  the 
data  are  plotted  as  axial  stress  versus  axial  strain,  the  slope 
of  which  Is  the  constrained  modulus  M.  When  the  data  are  plotted 
as  principal  stress  difference  versus  mean  normal  stress,  the 
slope  Is  2G/K  or.  In  terms  of  Poisson's  ratio  v,  3(l-2v)/l+v) . 

3)  The  third  type  of  UX  test  (designated  UX/Null)  Is  similar  to  the 
K0  test  In  that  both  radial  and  vertical  pressures  are 
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controlled.  A  wafer- shaped  specimen  Is  remolded  in  a  thin-walled  steel 
cylinder  which  Is  strain  gaged  on  the  outside.  As  vertical  pressure  is 
applied,  the  circumferential  strain  (measured  by  the  strain  gages)  on  the 
steel  cylinder  Is  kept  at  zero  by  applying  lateral  pressure  to  the 
-cylinder.  This  process  Is  continued  throughout  the  test.  The  data  are 
plotted  and  properties  deduced  the  same  as  those  for  the  UX/KQ  test. 
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APPENDIX  B 


Material  Model  Descriptions 

The  following  are  brief  descriptions  of  the  models  used  to  calculate 
,the  laboratory  behavior  of  the  sands  discussed  in  this  report. 

B.l  Elastic-Plastic  Model 

This  model  uses  a  linear  hydrostatic  stress-strain  relationship  in 
both  loading  and  unloading.  The  shear  stress-strain  response  is  also 
linear  up  to  a  failure  state,  after  which  the  material  behaves  as  an 
ideally  plastic  material.  The  failure  condition  is  defined  by  a  surface 
in  stress  space  which  is  symmetric  about  the  p-axis.  Figure  B-l  shows  the 
model  components  and  defines  the  model  parameters  as  used  in  the  SEM. 

B.2  AFWL  Engineering  Model 

This  model  is  essentially  a  version  of  the  elastic-plastic  model, 
modified  to  allow  piecewise  elastic  approximation  of  nonlinear  behavior 
and  volume  compaction.  The  model  is  defined  with  a  hydrostat  (in  pressure 
versus  volume  strain  space),  a  constant  Poisson's  Ratio  prior  to  failure, 
and  a  failure  surface  (in  invariant  stress  space).  Tensile  stresses  are 
limited  by  a  minimum  pressure  cutoff.  Plastic  strains  at  failure  are 
defined  by  a  nonassoclatlve  flow  rule,  because  stress  excursions  exceeding 
the  failure  surface  are  corrected  vertically  (l.e.,  at  constant  pressure) 
back  to  the  failure  state.  Figure  B-2  defines  the  SEM  parameters  for  the 
AFWL  model . 

B.3  Cap  Model 

The  behavior  of  the  cap  model  in  the  elastic  regime  is  governed  by 
elastic  bulk  and  shear  modulii.  The  bulk  modulus  is  a  function  of 
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pressure,  and  the  shear  modulus  depends  on  the  second  invariant  of  the 
stress  deviator  tensor,  The  plastic  behavior  of  the  model  is 

described  In  stress  space  by  a  modified  Drucker-Prager  failure  surface  and 
an  elliptical  strain-hardening  cap.  The  failure  surface  is  stationary  (in 
this  version),  and  the  movement  of  the  cap  is  controlled  by  the  amount  of 
plastic  volumetric  strain  the  material  has  experienced.  The  cap  model  is 
an  effective  stress  model.  In  that  it  can  represent  both  drained  and 
undrained  volumetric  behavior  simulating  effective  and  total  stress 
responses,  respectively.  Pore  pressure  response  of  the  material  is  then 
readily  calculated  by  subtracting  effective  from  total  stresses. 

Figure  B-3  Illustrates  the  cap  model  and  defines  the  SEM  cap  parameters. 
B.4  Hyperbolic  Model 

This  model  combines  a  hyperbolic  curve-fit  to  tri axial  shear  stress- 
strain  response  with  a  simple  bilinear  hydrostat  for  volumetric  response. 
Figure  B-4  shows  the  hydrostatic  and  shear  response  parameters,  and  how 
the  shear  parameters  are  calculated. 

B.5  Model  Parameters 

Tables  B-l  through  B-4  list  model  parameters  for  MB  and  RB  sands,  and 
RV  alluvium.  These  parameters  were  determined  on  the  basis  of  laboratory 
results  only. 
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Table  B-l.  Elastic-Plastic  Model  Parameters 


Parameter 

Name 

Units 

MB 

Sand 

(Dry) 

MB 

Sand 

(Sat) 

RB 

Sand 

(Sat) 

RB 

Sand 

(Sat) 

RV 

A1  luvium 
(0-6  m) 

9 

Bulk  Modulus 

BULK 

Pa 

156E6 

156E6 

387E6 

536E6 

188E6 

Shear  Modulus 

SHEAR 

Pa 

60E6 

60E6 

295E6 

206E6 

72E6 

CA 

Pa 

0.0 

0.0 

0.0 

50E6 

0.06E6 

Failure  Surface 

CB 

1/Pa 

0.0 

0.0 

0.0 

5.1E-9 

0.0 

CC 

Pa 

0.0 

0.0 

0.0 

50E6 

0.0 

CAM 

~ 

0.267 

0.257 

0.211 

0.0 

0.269 

Tensile  Strength 

TCUT 

Pa 

0.0 

0.0 

0.6E6 

0.0 

0.0 

Flow  Rule 

RULE 

- 

1.0 

1.0 

1.0 

1.0 

1.0 

Density 

RHOREF 

kg/m3 

1710 

2120 

1650 

2030 

1930 

Table  B-2.  Modified  AFWL  Engineering  Model  Parameters 


Parameter 

Name 

Units 

MB 

Sand" 

(Dry) 

MB 

Sand 

(Sat) 

RB 

Sand 

(Dry) 

RB 

Sand 

(Sat) 

RV 

Alluvium 
(0-6  m) 

Poisson's  Ratio 
(load) 

Poisson's  Ratio 
(uni oad) 

POISL 

0.35 

0.33 

0.35 

0.35 

0.33 

POISU 

0.25 

0.25 

0.25 

0.25 

0.28 

Loading 

BKO 

Pa 

85E6 

187E6 

175E6 

494E6 

209E6 

Bulk 

BKL1 

Pa 

191E6 

114E6 

518E6 

597E6 

153E6 

Modulus 

BKL2 

Pa 

346E6 

274E6 

439E6 

543E6 

330E6 

Unload  Bulk 

BKU 

Pa 

3600E6 

5850E6 

5622E6 

4537E6 

6667E6 

PI 

Pa 

2.2E6 

2.8E6 

1.7E6 

4.3E6 

2.3E6 

Breakpoints 

P2 

Pa 

14.8E6 

10.1E6 

11.6E6 

15.7E6 

20.0E6 

P3 

Pa 

30.0E6 

30.0E6 

30.0E6 

30.0E6 

53.0E6 

Tensile  Strength 

ST 

Pa 

0.0 

0.0 

-0.6E6 

0.0 

-0.1E6 

Intercept 

y 

Pa 

0.0 

0.0 

0.12E6 

0.0 

0.058E6 

FS  Slope 

s 

- 

0.267 

0.257 

0.211 

0.252 

0.269 

Von  Mlses  Limit 

VM 

Pa 

100E6 

100E6 

100E6 

50E6 

145E6 

Density 

RHOREF 

kg/m3 

1710 

2120 

1650 

2030 

1930 

Table  B-3.  Cap  Model  Parameters 


Parameter 

■» 

Name 

Units 

MB 

Sand 

(Dry) 

MB 

Sand 

(Sat) 

RV 

Alluvium 
(0-6  m) 

Effective 

AKI 

Pa 

208E6 

290E6 

160E6 

138E6 

200E6 

Stress 

AKI 

- 

0.978 

0.984 

0.969 

0.977 

0.983 

Hydrostat 

AK2 

1/Pa 

4.12E-9 

3.24E-9 

1.41E-8 

1.43E-8 

5.89E-8 

Total 

AKIM 

Pa 

- 

5840E6 

• 

8465E6 

5840E6 

Stress 

AKIM 

- 

- 

sm 

_ 

5.0E-4 

5.0E-4 

Hydrostat 

AK2M 

1/Pa 

- 

1.19E-8 

- 

7.22E-8 

1.20E-8 

Shear 

AGI 

Pa 

110E6 

170E6 

96E6 

65E6 

90E6 

Response 

AG1 

• 

0.900 

0.920 

0.900 

0.960 

0.990 

AG2 

1/Pa 

4.7E-6 

1.7E-7 

3.3E-6 

8.0E-6 

5.0E-6 

FS  Cohesion 

AC 

Pa 

0.0SE6 

0.0 

0.12E6 

0.1E6 

0.06E6 

FS  Slope 

AM 

• 

0.267 

0.257 

0.211 

0.252 

0.269 

Cap  Shape 

ARI 

- 

2.5 

3.0 

3.5 

3.5 

2.3 

Hardening 

AM 

- 

0.36 

0.36 

0.30 

0.30 

0.50 

AD 

1/Pa 

1.2E-8 

1.2E-8 

4.4E-9 

2.7E-9 

7.5E-9 

Density 

RHOREF 

kg/m3 

1710 

2120 

1650 

2030 

1930 

Table  B-4.  Hyperbolic  Model  Parameters 


Parameter 

Name 

Units 

MB 

Sand 

(Dry) 

Pfl 

Modulus  Number 

K 

990 

21 

Modulus  Exponent 

n 

- 

0.29 

1.45 

Failure  Ratio 

Rf 

- 

0.75 

0.99 

Cohesion 

c 

Pa 

0.0 

0.4E6 

Friction  Angle 

0 

degrees 

34 

26 

Unload/Reload  Modulus  No. 

Kur 

- 

24800 

10946 

G 

- 

0.44 

0.44 

Poisson's  Ratio 

F 

• 

0.08 

0.08 

d 

• 

8.2 

8.2 

Density 

RHOREF 

kg/m3 

1710 

1650 

APPENDIX  C 

Sample  Results  from  SEM  Finite  Difference  Calculations 


Notes 

Problem  No.  -  C22-4 

Geometry  -  One-D  Axi symmetric 

Material  -  RV  Alluvium 

Model  -  Laboratory  Based  cap 

Loading  -  CIST  (Exponential) 

Units  -  Meters,  Pascals,  Seconds 

TXX  -  Normal  Stress  (Pa) 

TYY  -  Hoop  Stress  (Pa) 

EX11  -  Impulse  (Pa-s) 

XVL  -  Velocity  (m/s) 

XDX  -  Displacement  (m) 

"Axial"  -  Radial  Direction 

"Radial"  -  Hoop  Direction 


Grid  Point 


Range  (m 


1 

0.30 

5 

0.91 

9 

1.52 

15 

2.44 

23 

3.66 

RRI)  SOIL  ELEMENT  MODEL  VERSION  1  .T1N03 


flRR  SOIL  ELEMENT  MODEL  VERSION  l  JRN83 


ORfl  SOIL  ELEMENT  MODEL  VERSION  t  JRN03 
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COMPARISON  OF  CONSTITUTIVE  MOOELS  SUBJECTED  TO  VARIOUS 
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SUMMARY.  Many  levels  of  constitutive  models  for  geological  materials  are  available  today  to  the 
practitioner,  from  very  complex  to  simple  linear-elastic.  As  a  result,  there  is  often  confusion  when  the 
time  comes  to  choose  a  model  which  will  give  the  appropriate  combination  of  simplicity  (and  therefore 
cost-effectiveness)  and  accuracy.  This  paper  describes  a  computer  code  which  has  been  developed  as  a  tool 
for  studying  material  constitutive  models.  The  code  allows  comparisons  of  model  behavior  and  parametric 
studies  to  determine  the  influence  of  model  omponents.  It  Is  Intended  that  the  "Soil  Element  Model",  as  it 
Is  called,  be  an  aid  In  model  development  and  also  in  choosing  a  mode)  which  best  suits  a  particular  boundary 
value  problem.  A  study  Is  presented  which  Illustrates  the  use  of  this  code  to  compare  the  ability  of  several 
material  models  to  replicate  laboratory  test  data.  The  test  data  (both  static  and  dynamic)  were  taken  in 
part  from  previous  efforts,  while  some  data  were  generated  specifically  for  this  purpose.  The  study 
concentrates  on  parametric  effects  and  on  Isolating  areas  for  model  Improvement  and  development.  In 
addition,  complex  loading  paths,  typical  of  those  which  are  Induced  by  blast  loadings,  are  imposed  on  the 
models  and  comparisons  are  made. 


INTRODUCTION 

There  Is  a  great  amount  of  computer  simulation 
currently  being  done  to  assess  the  response  cf 
soils  to  Impulsive  dynamic  loadings.  Finite 
difference  and  finite  element  calculations  are 
heavily  relied  upon  for  prediction  of  ground 
motions  and  soil-structure  Interaction  phenomena. 
Along  with  establishing  suitable  boundary 
conditions  and  the  mathematical  particulars  of  the 
finite  difference  (or  element)  code  Itself,  choice 
of  constitutive  relationships  and  parameters  seems 
to  be  most  influential  on  calculated  motions.  It 
Is  often  not  known  precisely  what  effect  choosing 
one  type  of  material  model  over  another  will  have, 
so  a  particular  type  of  model  is  usually  chosen 
based  on  past  experience  of  the  calculator. 
Specific  model  parameters  are  determined  by 
utilizing  as  much  laboratory  and/or  In  situ  data 
as  is  necessary  and  available.  The  process  of 
fitting  parameters  to  data  may  range  from  quick 
and  easy  to  tedious  and  time  consuming,  depending 
on  the  complexity  of  the  model  and  the  amount  of 
experience  the  "user"  has  had  with  the  model. 

In  an  attempt  to  address  some  of  the  questions 
regarding  dynamic  soil  modeling,  a  computer  code 
called  the  Soil  Element  Model  has  been  developed 
(X).  The  prime  motivation  behind  the  Soil  Element 
Model  is  to  allow  comparison  of  different 
constitutive  models  when  they  are  subjected  to 
various  kinds  of  statle  and  dynamic  stress  and/or 


strain  paths.  In  addition,  the  code  allows  the 
user  to  become  familiar  with  the  response  of  each 
installed  model  through  slaple  parameter  variation 
studies.  Finally,  the  Soil  Element  Model  can  be 
utilized  as  a  model  fitting  tool.  It  facilitates 
the  iterative  fitting  of  parameters  to  data  which 
Is  required  by  some  models  and  can  be  used  to 
determine  the  sensitivity  of  model  response  to 
small  parameter  variations. 

SOIL  ELEMENT  H00EL 

Description 

The  Soil  Element  Model  (SEM)  consists  of  a 
driver,  a  set  of  boundary  condition  algorithms, 
and  a  set  of  constitutive  models.  Fig.  1 
Illustrates  the  basic  logic  of  the  code. 

The  program  Is  written  in  FORTRAN  and  is 
adaptable  to  most  any  computing  system.  The 
examples  for  this  paper  were  done  with  a  Cray  1-S 
computer  at  the  Air  Force  Weapons  Lao  (AFWL),  but 
the  code  Is  also  functioning  on  an  in-house  Apple 
II  microcomputer. 

Capabilities 

Three  types  of  boundary  conditions 
representing  both  laboratory  tests  and  in  situ 
conditions  may  be  applied:  stress  controlled, 
strain  controlled,  and  the  equations  of  motion 
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Fig.  1.— Basic  Soil  Element  Model  Logic 


governing  one-dimensional  wave  propagation.  At 
present,  the  stress  boundary  options  arc 
hydrostatic  compression,  proportional  loading,  and 
standard  trl axial  shear.  The  strain  boundary 
options  are  uniaxial  strain,  pure  shear,  true 
trlaxlal  (rigid  platen),  and  arbitrary  strain 
paths. 

Strain  controlled  boundary  conditions-  are 
easily  enforced  because  all  of  the  models  are 
formulated  such  that  given  a  strain  Increment, 
they  calculate  a  stress  Increment.  Stress 
boundaries  require  an  Iterative  scheme  to 
establish  a  suitable  strain  Increment  which  will 
produce  the  correct  stnmss  Increment.  When 
simultaneous  Iteration  Is  necessary  for  more  than 
one  strain  increment,  the  process  becomes  quite 
time  consuming.  (A  true  trlaxlal  test  with 
flexible  platens,  for  example,  may  require 
determination  of  three  independent  strain 
increments. ) 


Exercising  a  material  model  over  the  above 
stress/straln  paths  provides  Insight  Into  model 
response  under  many  conditions.  However,  since 
the  ultlMte  use  of  these  models  is  In  dynamic 
calculations,  a  relatively  simple  one-dimensional 
wave  propagation  code  was  Incorporated  Into  the 
SIN  at  a  boundary  condition  option.  This  allows 
the  user  to  establish,  at  least  In  one-dimension, 
whetner  differences  observed  in  lab  tests  or 
(Q^othetlcal  In  situ  model  behavior  produce 


significant  differences  in  an  actual  calculation. 
The  explicit  finite-difference  scheme  which  was 
lavlamented  Is  an  adaptation  of  a  more  general 
code  called  SNEAKY  written  by  Hart  (2). 
Calculating  wave  propagation  enforces  conservation 
of  mass  and  momentum  between  stacked  elements  (or 
zones).  It  requires  maintaining  a  separate 
constitutive  relation  reflecting  individual  stress 
and  strain  histories  (for  each  zone). 

Constitutive  Models 

Although  the  SEM  is  used  primarily  to  study 
soil  modeling,  the  code  Is  designed  to  accept 
arbitrary  constitutive  models  which  may  be  used 
for  any  material,  not  just  soil.  In  fact,  many  of 
the  models  studied  In  conjunction  with  dynamic 
soil  behavior  have  been  either  developed  for  or 
used  to  model  other  materials  such  as  concrete  and 
rock. 

Models  which  are  presently  Incorporated  In  the 
SEM  Include  elastic,  viscoelastic,  elastlc- 
perfectly  plastic,  one-dimensional  curve-fits 
(3,4),  AFVA  Engineering  (5),  and  the  Cap  model  (6). 

The  elastic  constitutive  relationship  Is 
simply  Hooke's  Law  for  a  linearly  elastic. 
Isotropic  materia).  The  model  contains  two 
constants,  a  bulk  modulus  (K)  ana  a  sncar  modulus 
(6),  although  any  other  two  elastic  constants  may 
be  substituted. 
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The  three-element  viscoelastic  model  consists 
of  •  spring  (Kj )  in  parallel  with  a  dashpot 
(nl),  both  In  sarlas  with  a  second  spring 
(Kj).  The  governing  equation  for  this  modal  may 
be  used  to  describe  either  volumetric  or  shearing 
behavior  for  soil.  It  may  also  be  Included  In 
other  models  to  account,  when  necessary,  for  rate 
dependent  elastic  behavior 

For  the  elastlc-perfectly  plastic  model,  a 
modified  Drucker-Prager  (pressure  dependent) 
failure  surface  Is  defined  which  limits  the 
devlatorlc  component  of  stress  and  allows  plastic 
strain  according  to  an  associated  flow  rule.  For 
stress  states  below  the  failure  surface  the 
material  Is  linearly  elastic.  In  all,  six 
parameters  are  required,  two  for  elastic  behavior 
and  four  for  the  exponentially  curved  failure 
surface  with  tensile  cutoff. 

One-dimensional  curve-fitting  Is  a  comnon 
empirical  method  for  modeling  soil  behavior. 
Pyke’s  curve-fit  (3)  for  shearing  behavior  fits  a 
hyperbolic  expression  to  the  initial  loading  shear 
stress-shear  strain  curve,  then  uses  extended 
Maslng  rules  to  model  the  hysteresis  loops.  At 
present,  a  simple  linear  elastic  relation  Is  used 
In  the  SEM  to  describe  the  volumetric  response, 
but  It  Is  possible  to  combine  Pyfcc's  shear 
behavior  with  more  complicated  bulk  responses  for 
better  overall  behavior. 

The  hyperbolic  curve-fitting  method  of  Ouncan 
and  Chang  (4)  collapses  trlaxlal  shear 
stress-strain  results  quite  welt.  In  conjunction 
with  a  confining  pressure-dependent  bulk  modulus, 
the  hyperbolic  relations  can  be  used  for  strain 
predictions.  Because  the  relations  are  direct 
representations  of  data,  the  lab  response  of  many 
different  types  of  soil  may  be  adequately 
modeled.  The  models  lack  the  capability  to  match 
more  complex  stress/straln  paths,  which  may  or  may 
not  be  critical,  depending  on  the  particular 
application  of  the  hyperbolic  models. 

The  AFWL  Engineering  model  may  be  classified 
as  elastlc-perfectly  plastic  with  volume 
hysteresis.  Behavior  Is  defined  by  a  plecewlse- 
Unear  hydrostat  and  a  Drucker-Prager  failure 
surface  with  a  von-HI ses  cutoff.  The  model  Is 
commonly  used  In  ground  shock  calculations.  The 
SEM  version  has  twelve  parameters  wnlch  are 
relatively  easy  to  determine  from  standard 
laboratory  tests  (uniaxial  strain  and  trlaxlal 
shear),  or  they  may  be  estimated  from  In  situ  test 
data. 

The  Cap  model  Is  also  commonly  used  for 
calculating  ground  shock  effects.  It  is  available 
on  several  different  levels  of  complexity  which 
account  specifically  for  various  behaviors,  such 
as  rate  dependence,  anisotropy,  kinematic 
hardening,  etc.  Anyone  wishing  to  use  the  Cap 
model  must  first  decide  which  behavioral  phenomena 
are  Important  to  their  particular  problem  and  then 
be  sure  that  the  correct  features  are 
Implemented.  The  version  currently  Installed  in 
the  SEM  Is  an  effective  stress  model  with 
llneai^elastlc  behavior  within  the  failure  surface 
and  cap,  a  curved  or  modified  Orucker-Prager 
failure  surface,  and  an  elliptical  cap.  The 
failure  surface  is  fixed  and  movement  of  the  cap 
is  controlled  by  plastic  volumetric  strain.  To 
define  both  drained  (total  stress)  and  undrained 
(affective  stress)  behavior,  16  parameters  must  be 


determined.  A  uniaxial  strain  or  hydrostatic 
compression  test  and  both  drained  and  undrained 
trlaxlal  shear  tests  (with  pore  pressure 
measurements)  are  necessary. 

Additional  models  which  can  accomodate  more 
complex  soil  behaviors  are  currently  being 
installed  in  the  SEM.  Composite  models,  which 
combine  various  desirable  features  of  different 
models,  are  also  oeing  studied. 

SEM  APPLICATION 

The  following  section  is  an  example  of  how  the 
Soil  Element  Method  is  being  applied  to  better 
understand  and  improve  the  modeling  of  dynamic 
soil  behavior. 

Material  Description 


The  SEM  exercises  described  here  are  aimed  at 
evaluating  the  response  of  Misers  Bluff  (MB)  sand, 
a  medium  to  coarse  grained  silty  sand  (Sw-SM) 
obtained  at  Planet  Ranch,  Arizona.  Planet  Ranch 
was  the  site  for  a  cylindrical  in  situ  test  in 
197?  (CIST  19;  8)  and  a  series  of  high  explosive 
events  In  1978  (Misers  Bluff  II;  9).  For  these 
modeling  studies,  available  lab  data  for  the 
behavior  of  saturated  MB  sand  was  collected  and 
additional  tests  were  performed  on  air  dry  samples 
remolded  to  a  density  of  1.72  g/cc.  The  net 
result  was  a  fairly  complete  set  of  data  for  the 
uniaxial  strain  and  triaxlal  shear  responses  of 
wet  and  dry  MB  sand  to  stress  levels  as  high  as  60 
MPa  and  confining  pressures  up  to  6.9  MPa  (see 
Phillips;  10). 

Laboratory  Behavior 

Parameters  for  the  elastic-plastic,  AFWL 
Engineering,  and  Cap  models  were  obtained  by  using 
the  SEM  to  fit  the  model  responses  to  the 
corresponding  representative  laboratory  data.  The 
elastic-plastic  and  AFHL  models  Involve  straight¬ 
forward  linearization  of  data,  while  the  Cap 
requires  some  Iteration  to  find  the  right 
combination  of  parameters  to  match  both  the  bulk 
and  shear  behaviors. 

As  shown  In  Fig.  2,  the  AFWL  and  Cap  models 
can  match  the  nonlinear  uniaxial  strain  response 
fairly  well,  while  the  straight  line  approximation 
of  the  elastic-plastic  1$  probably  not  acceptable 
for  most  purposes.  At  stress  levels  below 
failure,  the  elastic  model  is  not  able  to  produce 
irrecoverable  deformations.  The  Cap  model  is  the 
most  capable  of  the  three  for  matching  the 
laboratory  shear  behavior,  as  seen  In  Fig.  3. 
Note  that  the  Cap  model's  trlaxlal  stress-strain 
response  could  be  brought  even  closer  to  the  data 
with  more  Iteration. 

All  three  models  have  a  single  pressure 
dependent  failure  surface  with  a  circular 
Intersection  in  the  octahedral  plane.  This  does 
not  compare  favorably  with  the  data  on  MB  send, 
which  Indicates  asymmetry  of  the  failure  surface 
about  the  P-axIs.  A  Nohr-Coulomb  shape  would 
probably  be  more  representative.  An  associative 
flow  rule  in  the  elastic-plastic  and  Cap  models 
allows  shear-dl latency  coupling  to  take  place, 
while  the  AFWL  model  uses  a  nonassoclatlve  flow 
rule  on  the  failure  surface  and  has  no  coupling. 
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In  Situ  Behavior 

It  It  util  known  that  toll  behavior  observed 
In  the  flald  cannot  alwayt  be  well  predicted  by 
constitutive  aodels  bated  on  laboratory  tatt 
ratultt.  The  ttratt/ttraln  patbt  encountered 
In  tltu  are  often  aucb  different  and  sore  complex 
than  those  In  atandard  laboratory  tattt, 
atpeclally  In  the  cate  of  lapulslve  loadt.  In 
addition,  there  are  the  probloat  Involved  with 
sampling  sol It  which  may  create  significant 
differences  in  lab  versus  in  situ  properties. 
There  are  various  ways  of  coping  with  these 
proolaas.  first,  one  could  assuaa  the  lab  Model 
Is  adequate  and  proceed  with  calculating.  Second, 
the  laboratory-derived  paraaeters  could  be 
adjusted  to  better  reflect  In  situ  response. 
Adjustments  aay  be  based  on  experience  and/or 
Intuition.  Lab  tests  could  be  performed  which  try 
to  slaulate  In  situ  stress  paths  aore  closely 
(e.g.  true  trlaxlal).  Or  finally.  In  situ  data 
aay  be  taken  and  lap  roved  node  Is  could  be 
formulated  based  on  these  results,  with  each  of 
these  alternatives  coaes  increased  expense  and, 
usually,  aore  complex  constitutive  relationships. 

A  recent  example  of  an  atteapt  to  slapllfy  the 
In  situ  aodellng  process  Is  the  hypothesis  of 
Morkaan,  at  al  (11)  that  strain  paths  In 
explosively  driven  fields  of  notion  are  Halted  to 
a  snail  region  In  strain  space.  These  dominant 
strain  paths  are  proposed  to  be  insensitive  to 
variations  In  geology  and  primarily  dependent  on 
shot  geometry.  If  material  behavior  could  be 
measured  in  the  lab  over  these  paths  (Indications 
are  that  this  would  require  servo  strain- 
controlled  apparatus),  modeling  would  be  greatly 
simplified.  In  fact,  the  measured  stress  path 
would  effectively  become  the  constitutive  model. 


Fig.  3.— Misers  Bluff  Sand  Trlaxial  Shear  Response 


Fig.  4  Is  one  hypothetical  stress  path  for  a 
spherically  symmetric  ground  motion  field.  Note 
that  the  strain  Is  nearly  uniaxial  until  the  point 
where  radial  relief  occurs.  This  strain  path  was 
used  In  the  SEM  to  drive  the  elastic-plastic, 
AFVC,  and  Cap  models  for  MB  sand.  Three  levels  of 
Initial  K0  consolidation  were  applied  to  study 
how  predicted  behavior  might  vary  with  depth. 
Results  for  the  three  models  and  three  confining 
pressures  are  compared  In  Fig.  5.  The  differences 
In  resultant  stress  paths  reflect  to  some  extent 
the  same  differences  seen  In  lab  shear  behavior 
but  also  reflect  basic  model  differences.  Shear 
failure  occurs  upon  unloading  due  to  radial  relief 
and  the  failure  surface  influences  all  cases. 
With  depth,  the  K0  condition  diverges  from  the 
failure  surface  and  the  material  undergoes 
Increased  stress  relief.  (Tensile  relations  are 
also  quite  Important  and  effect  the  late  time 
predicted  response. ) 

The  relative  differences  In  predicted  stress 
paths  may  Indicate  that  strain  paths  saen  in  a 
calculation  are  In  fact  dependent  on  choice  of 
constitutive  relationships  and  may  also  vary  with 
geology  as  determined  by  specific  material 
parameters.  Further  SEM  studies  will  be  helpful 
In  assessing  the  dominant  strain  path  hypothesis. 

In  situ  models  are  sometimes  determined  by 
Iteratively  performing  finite  difference 
calculations  while  varying  model  parameters  until 
the  data  waveforms  are  matched.  The  Inclusion  of 
SNEAKY  in  the  SEM  facilitates  the  evaluation  of 
various  models  in  this  context.  Fig.  6  shows  the 
setup  for  a  one-dimensional,  one  layer  simulation 
of  a  DISC  TEST  event  (12).  The  MB  sand  models 
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Fig.  Characteristic  Strain  Path  for  Spherically 
Symmetric  Ground  Notion  (Ref.  4,  p.  70) 


used  In  the  previous  examples  were  again  used. 
Simulation  time  Is  short  (40  ms)  and  loading  1$ 
typical  of  data  from  the  DISC  TEST  I  and  II  events. 

Velocity  waveform  comparisons  are  shown  In 
Fig.  7a.  For  the  most  part,  differences  In  peak 
velocities  and  waveform  shapes  are  due  to 
differences  In  uniaxial  strain  behavior  of  the 
models  below  10  IPa  (Fig.  2).  Note  the  similarity 
of  the  AFWL  and  Cap  waveforms,  the  two  models 
perform  similarly  in  one-dimensional  wave 
propagation.  Attenuation  of  peak  velocity  with 
depth  (Fig.  7b)  Is  slightly  faster  for  the  AFWL 
model  than  for  the  Cap.  Evidently  the  curved 
nature  of  the  Cap  model's  unload  volumetric 
relation  has  an  appreciable  effect  on  energy 
dissipation.  As  expected,  the  elastlc-plstlc 
model  exhibits  almost  no  attenuation  with  depth. 

CONCLUSIONS 

The  complexity  of  each  calculatlonal  problem 
dictates  the  nature  of  the  constitutive 
relationship  which  must  be  used  to  accurately 
model  material  response.  The  available  models  are 
eacn  adequate  within  the  limits  of  the  behaviors 
they  were  formulated  to  deal  with.  When  these 
bounds  are  exceeded,  the  results  from  any  model 
should  be  viewed  with  skepticism.  An  Idea)  model 
from  the  practitioner’s  viewpoint  would  be  capable 
of  reproducing  all  soil  behaviors,  but  would  be 
formulated  such  that  unnecessary  effects  could  be 
switched  off.  In  this  way,  the  modeler  would  be 
using  the  simplest  model  possible  for  each 
particular  problem  while  maintaining  consistency 
from  one  application  to  another. 


Fig.  5.— Resultant  Stress  Paths 


The  SEN  Is  Intended  to  be  used  as  a  tool  for 
evaluating  the  bounds  of  a  constitutive 
relationship,  for  both  laboratory  ana  more 
complicated  In  situ  cases.  This  approach  is 
necessary  to  make  rational  choices  when 
considering  the  application  of  any  particular 
constitutive  model. 
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Fig.  7.— Comparison  of  Computed  Motions: 
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with  Depth 


(8)  Amend,  J.H.,  and  Brown,  R.W.,  "Cylindrical 
In  Situ  Tests  at  Selected  Test  Sites",  Air 
Force  Weapons  Laboratory,  Klrtland  Air  Force 
Base,  NM,  August,  1980. 

(9)  Strode,  J.D.  Jr.,  et  al,  "Proceedings  of  the 
Misers  Bluff  Phase  II  Results  Symposium, 
Vols.  I-III",  Report  No.  POR  7013-1,  Field 
Coanand,  Defense  Nuclear  Agency,  Klrtland 
Air  Force  Base,  NM,  March  1979. 

(10)  Phillips,  B.R.,  "Mechanical  Properties  of 
Misers  Bluff  Sand",  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS, 
May,  1982. 

(11)  Workman,  J.W.,  Trulio,  J.G.,  and  Stokes, 
E.S.,  "Modeling  the  Behavior  of  Geologic 
Materials  in  Explosive  Field  Events",  Report 
No.  AFWL-TR-80-66,  Applied  Theory,  Inc., 
Los  Angeles,  CA,  January,  1981. 

(12)  Jackson,  A.E.  Jr.,  Murrell,  D.W.,  Zelasko, 

J.S.,  and  Skinner,  F.W.  Jr.,  "Ralston  valley 
Soil  Compressibility  Study:  Quick-Look 

Report  for  OISC  TEST  I",  U.S.  Army  Engineer 
Waterways  Experiment  Station,  Vicksburg,  MS, 
May,  1981. 


SOIL  MODEL  EVALUATION  UNDER  DYNAMIC  LOADINGS 


by 

William  C.  Dass1 
Jimmie  L.  Bratton2 


Abstract 

There  are  many  different  types  of  constitutive  relationships 
available  for  calculating  the  response  of  geologic  media  to 
impulsive  dynamic  loading.  Choosing  a  material  model  which  is 
suitable  for  a  particular  situation  can  be  difficult.  Model 
selection  Is  often  based  on  personal  experience  of  the  calculator 
which  may  or  may  not  result  In  optimum  efficiency  and  best  response. 

This  paper  describes  a  computer  code  which  has  been  developed  as  an 
aid  for  studying  material  constitutive  models.  The  Soil  Element  Model 
can  calculate  the  response  of  a  given  material  model  to  laboratory 
and  In  situ  test  conditions,  arbitrary  strain  paths,  or  one-dimensional 
wave  propagation.  It  Is  useful  for  developing  models,  performing  parametric 
studies  to  determine  model  component  Influence,  and  comparing  model  behaviors. 
A  study  Is  presented  which  Illustrates  the  use  of  this  code  to  compare  the 
ability  of  several  material  models  to  replicate  laboratory  and  In  situ 
data.  The  study  focuses  on  a  sand  from  the  DRY  CARES  site  near  Yuma, 

Arizona,  and  examines  the  advantages  and  disadvantages  of  each  model  selected. 
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